

1 -60-0184 697 DETERMINATION 

1 CONJUGATION IN 

SCHOOL NONTERE 

| UNCLASSIFIED 

OF THE 
BARIU 
* CA 

ORIGI 
H TITfi 
T R MG 

N OF SELF-PUMPED PHASE 
NATECU) NAVAL POSTGRADUATE 
ORE JUN 87 

F/G 26/2 

1/ 

Nl 


i K 

■ 



























n 































■ 

■ 

" 








- 




"n 

m 

m 

S3 

SI 

a 


m 

m 

■ 


Rl 
















_ " --— ___:_ J 






























AD-A184 697 


NAVAL 


POSTGRADUATE SCHOOL 

Monterey, California F|[[ 





THESIS 


DETERMINATION OF THE ORIGIN OF 
SELF-PUMPED PHASE CONJUGATION IN 
BARIUM TITANATE 

by 

Thomas R. Moore 

June 1987 

Thesis Advisor Donaid L. Wallers 

Approved for public release; distribution is unlimited. 


77 7 <33 33 i 









lj REPORT SECURITY CLASSIFICATION 

UNCLASSIFIED 


2a security Classification authority 


REPORT DOCUMENTATION PAGE 


ib. restrictive markings 


2b DECLASSIFICATION/DOWNGRADING SCHEDULE 


4 PERFORMING ORGANIZATION REPORT NUMBER(S) 


J DISTRIBUTION'/AVAILABILITY OF REPORT 


Approved for public release; 
distribution 


S MONITORING ORGANIZATION REPORT NUMBER(S) 


6a NAME OF PERFORMING ORGANIZATION 6b OFFICE SYMBOL 7a NAME OF MONITORING ORGANIZATION 

(If applicable) 

Naval Postgraduate School 61 Naval Postgraduate School 


6c AOORESS (C/ry. Stare, and 21* Cod*) I 7b ADDRESS (City, Star*..and II* Cod*) 


Monterey, California 93943-5000 


8a NAME OF FUNOiNG/SPONSORING 
ORGANIZATION 


Monterey, California 93943-5000 


8b OFFICE SYMBOL I 9 PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 
(if ippHabl*) j 


10 SOURCE OF FUNDING NUMBERS 


PROGRAM 
ELEMENT NO 


PROJECT - 

TASK 

NO 

NO 


8c AOORESS (Cry. Set*. and it*Cod*) 


H TiTlE (ln<lud* Security Clsuificstion) 

DETERMINATION OF THE ORIGIN OF SELF-PUMPED PHASE CONJUGATION IN BARIUM TITANATE 


U personal auTmoR(S) 


WORK UNIT 
ACCESSION NO 





lib TIME COVEREO 
FROM _ TO 


14 OATE OF REPORT (Ytsr, Month. Diy) IS PAGE COUNT 

1987 June 98 



COSATi COOES 


GROUP SUBGROUP 


18 SUBJECT TERMS (Continue on reverie if neceuary and identify by biock numbtr) 
Phase-conjugation, Barium Titanate, Photorefractors 
Stimulated Photorefractive Scattering 


•9' ABSTRACT (Continue on reverie if neceuary and identify by bloc* numbtr) 

\\ 'Bar lUHOh ,rtC' 

When a beam of extraordinarily polarized light is directed through a 
poled, single crystal of CjJaTlO^it fans in the direction of the optic 
axis. This thesis research "investigates the origin of this fanning, or 
asymmetric self-defocusing, in barium titanate. Experimental 
observations and theoretical simulation suggest that this phenomenon is 
due to near forward stimulated ^photorefractive s scattering. 
Additionally, observations of unique variegated beam patterns during 
phase conjugation are reported, including internally pumped rings and 
motion invariant patterns. 


20 0 $TRi9UTiON/AVAILA8ILiTY OF ABSTRACT |\ 21 ABSTRACT SECURITY CLASSIFICATION 

BuNCLASSlFlEOAJNLiMlTEO □ SAME AS RPT ! □ OTIC USERS JNCLASSIFIED 


22a NAME OF RESPONSIBLE 'NOlVlOUAL 

T H T.. TJa 1 


DO FORM 1473, 84 MAR 83 APR edition may be used unt.l exhausted SECURITY CLASSi 

All other tditiont art obtolat# 


SECURITY CLASSIFICATION OF This PAGE 












































Approved for public release; distribution is unlimi ted. 


Determination of the Origin of 
Self-Pumped Phase Conjugation in 
Barium Titanate 


by 


Thomas R. Moore 
Captain, United States Army 
B.S., Stetson University, 1978 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN PHYSICS 


from the 


NAVAL POSTGRADUATE SCHOOL 
June 1987 



2 







ABSTRACT 


When a beam of extraordinarily polarized light is directed through a poled, single 
crystal of BaTi0 3 it fans in the direction of the optic axis. This thesis research 
investigates the origin of this fanning, or asymmetric self-defocusing, in barium 
titanate. Experimental observations and theoretical simulation suggest that this 
phenomenon is due to near forward stimulated photorefractive scattering. 
Additionally, observations of unique variegated beam patterns during phase 
conjugation are reported, including internally pumped rings and motion invariant 
patterns. 
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I. INTRODUCTION 


During the last decade the photorefractive crystal barium titanate (BaTiO,) has 

L---- , 3 

received a large amount of attention within the scientific community. Barium titanate 
has long been of interest to acousticians and solid state physicists because of it's 
photoacoustic and piezoelectric properties, but only recently have optical physicists 
become very interested in this crystal. The cause of this sudden and intense interest 
was the discovery of photorefraction and optical phase conjugation. Originally found 
in liquids and gases [Ref. 1], optical phase conjugation has the ability to reconstruct a 
distorted wavefront (sometimes loosely referred to as time reversal). Barium titanate is 
one of many photorefractive crystals which can exhibit phase conjugation. 

There exist a number of intriguing optical properties of single crystals of barium 
titanate, but by far the most interesting is its ability to phase-conjugate at continuous 
wave low powers, with no external apparatus. This effect, known as self-pumped phase 
conjugation (SPPC), is the subject of extensive current investigation. The goal of this 
research is to determine the origin of self pumped phase conjugation in BaTi0 3 . 

Currently there exist two explanations of the cause of SPPC in BaTiOj. One 
theory, developed by Feinberg and his students at UCLA [Ref. 2], attributes the 
production of the phase conjugate reflection to the photorefractive effect via 
Degenerative Four Wave Mixing (DFWM). According to this theory, an asymmetric 
change in the index of refraction across the path of an incident beam causes the beam 
to bend. This bending eventually leads to a geometry which, through two internal 
reflections (often termed a CAT comer), allows the beam to intersect itself, leading to 
DFWM. The second theory, originating with Lam and his co-workers at Hughes 
Research Laboratory [Ref. 3], attributes SPPC to Simulated Photorefractive Scattering 
(SPS), or Two Wave Mixing, similar to that exhibited in gases and liquids with two 
frequency disparate input beams of very high powers. With the goal of proving one or 
the other of these theories, this thesis begins with a general background and theory of 
phase conjugation as it applies to BaTi0 3 (Chapter II). Chapter III describes PSIM 
(Photorefractive SIMulation), a simulation of the phenomenon of beam fanning due to 
a change in index of refraction within the incident beam path arising from the 
photorefractive effect. This simulation enables the prediction of the change in index of 







refraction and the index profile necessary to account for SPPC with only the 
photorefractive effect taken into account. 

Chapter IV introduces experiments to determine the effect of incident beam size 
on various crystals of BaTiOj. Chapter V summarizes the results and analyzes the 
data presented. The conclusion is reached that Stimulated Photorefractive Scattering :s 
responsible for the process of beam fanning which, leads to self pumped phase 
conjugation in BaTiOj. The actual origin of the phase conjugate beam is attributed to 
both SPS and DFWM. 

Additionally, Appendix A contains a summary of the work I did at Los Alamos 
National Laboratory under the direction of Dr. R. A. Fisher in February and March, 
1986. This work describes and catalogs internal beam production in barium titanate 
during SPPC. Some internal beam patterns not previously reported are introduced and 
explanations for these phenomena are proposed. The results of the experimental work 
in the body of this thesis have been submitted for publication to Optics Letters. An 
expanded version of Appendix A was presented at the XIV International Quantum 
Electronics Conference (1986), and has been submitted for publication in the Journal of 
the Optical Society of America B in collaboration with R. A Fisher and A. V. Nowak. 




II. BACKGROUND AND THEORY 


A. INTRODUCTION 

Optical phase conjugation, first discovered in 1972 [Ref. 1], relies on a specific 
nonlinear property of a substance to create the formation of the complex conjugate of 
an incident wave. The reflection involved is not the familiar specular reflection seen in 
everyday life, but a reflection that retraces the path of the incident beam exactly, 
regardless of the direction or angle of incidence. The exact retracing of the incident 
beam path and the formation of the complex conjugate of the incident wavefront are 
the two unique and most important effects that make phase conjugation both 
interesting and useful. In barium titanate this phase conjugate beam is directly a result 
of the photorefractive effect. 

B. ORIGINS OF THE PHASE-CONJUGATE BEAM 

The origin of the phase conjugate beam in photorefractors is light which is Bragg 
scattered by an index grating formed by the light in the material. This process is 
known as wave mixing. There exists both two-wave mixing (scattering) and four-wave 
mixing in barium titanate. Both of these processes will be treated theoretically in 
general terms; the specifics of wave mixing in BaTiOj will be treated later. 

To understand the origins of phase conjugation one must begin with the wave 
equation for a homogeneous, isotropic, time-independent, nonmagnetic and 
nonconducting dielectric. The polarization term is divided into its linear (Pj) and 
nonlinear (P n) ) pans, and the linear portion is absorbed into the linear permittivity of 
the medium. £. This leaves the wave equation as: 

, d 2 E d 2 P„, 

V 2 E -£H— -p—f 41 . (2.1) 

ct" Sv 

For ease of manipulation an approximation to this equation is made that reduces 
the second order nonlinear wave equation to a first order equation. 1 This 
approximation is derived by Fisher and Yariv [Ref. 5: pp. 9-11] and for it to be valid, 
the envelope of the pulse must not change appreciably during an optical period, hence 

l Some calculations can be made without this approximation but they are quite 
tedious (e.g. [Ref. 4]). 







the term Slowly Varying Envelope Approximation. In making the approximation it is 
necessary to assume that the electromagnetic waves under consideration are plane 
waves, and also to extract only the portion of the nonlinear polarization that 
propagates with the same frequency and wave vector as the modulation function of the 
electric field. This, referred to as the phase matched portion of the polarization, is 
given by 

Pni s P( z >t) ex P^i(tot ” kz), (2.2) 

in the wave given by 


E s E(z,t) exp ± i(cot — kz). (2.3) 

Here the boldface letters (E and P) are the fields and the conventional letters (E and P) 
are the field envelopes, which vary little in an optical period and in an optical 
wavelength. All other terms are ignored because they are not phase matched and 
therefore cannot couple to the electric field. The Slowly Varying Envelope 
Approximation, or SVEA, is then invoked. That is: 

, dE d 2 E 

|k 2 E| >>|k—1> >1-2-1. (2.4) 

oz oz i 

Fisher and Yariv then show that the wave equation with SVEA becomes a first-order 
equation relating the envelope functions E and P: 


dE . dE 
—- + Jz\i— 

OZ Cl 


itoVM 

±— ~ P. 

2 Vs 


(2.5) 


This result is the SVEA wave equation for a plane wave traveling in the - z direction, 
which snows that, given a nonlinear polarization, an eiectnc field can be produced. 

This development can be generalized to the case of many waves present within a single 
medium, resulting in the origin of the nonlinear processes that can lead to phase 
conjugation. 

The polarization term can often be written as a power series expansion: 


12 





P - EX{E) - X^E + X (2) E 2 + X (3) E 3 


( 2 . 6 ) 


where X(E) is the susceptibility of the material in which the wave is traveling. The 
expansion of the susceptibility in a power series separates the effects of the nonlinear 
material into distinct groups that categorize nonlinear effects [Ref. 5: p. 13j: 

• X^ - These are the linear properties that are the subject of classical optics and 
are not directly responsible for nonlinear effects. 

• X^ * The second order effects are in general called three-wave mixing and 

include: second-harmonic generation, optical rectification, parametric mixing 
and the Pockels effect. 

• X^ 3 * - Third order effects include some of the most popular effects responsible 
for phase conjugation such as: third-harmonic generation, nondegenerate four- 
wave mixing, Raman scattering, ac Kerr effect (degenerate four-wave mixing), 
dc Kerr effect, Brillouin scattering, two photon absorption, dc-induced 
harmonic generation and Stimulated Photorefractive Scattering. 

The two methods of creating a phase-conjugate reflection that are directly applicable to 

the experimentation described here are scattering and degenerate four wave mixing. 

1. Four-Wave Mixing 

In 1977 Hellwarth [Ref. 6] showed that it was possible to generate a time- 
reversed replica of any monochromatic-beam wave pattern. The method employed the 
interaction of an incident beam with counter-propagatmg pump waves in a 
homogeneous, transparent, nonlinear medium (see Figure 2.1). This method, called 
four-wave mixing, has become one of the most popular methods for phase conjugation 
because of the almost universal availability of the necessary equipment and proper 
media. 

The description of wave mixing begins with the nonlinear polarization derived 
from the energy function. As shown by Pepper and Yariv [Ref. 5: pp. 26-33). isolating 
only the third order term gives 



4\ 3 F F F 




. vOnsiuer acre regenerate lo^r-.vu'.e mixing vDl 'A Ni , a nere 
(Dj =■* = CD? = of represents frequencies of the four beams The subscripts 1 and 2 
refer to two strong, precisely counterpropagatmg pump beams, the subscript p refers to 
the weak probe beam, and the subscript c refers to the phase-conjugate output wave. 
It is assumed that the following conditions are met: 









Figure 2.1 Schematic of the process of four-wave mixing [Ref. 5: p.4'21], 

a) intense (nondcpletcd) counterpropagating planewave pump beams of equal 
intensity (kj, k,) t 

b) a weak probe beam (k p ), 

c) both the intense pump beams and the weak probe beam are at the same 
frequency to. 

Since the two pump waves are very strong compared to the probe beam they form a 
standing wave. To obtain the steady-state solutions, we can set 



( 2 . 8 ) 


oo. assuming ;ae $VE.\ and considering only the phase matched portion of the phase 

conjugate wave, then from Equation 2.5 












Where 


“ 77 x ’ e > e 2 


(2.16) 


E C (L) is the amplitude of the conjugate wave at its origin (the phase conjugate mirror), 
and E p (0) is the amplitude of the probe wave at its origin. 

In practice, the amplitude of the phase conjugate wave at its origin is zero (i.e. 
E c (L) = 0), so that at the input point (z= 0) of the probe beam 

k * 

E c (0) - i — tan(|k|l)E p *(0). (2.17) 


Equation 2.17 shows that at z = 0 (the origin of the probe beam) the reflected field is 
proportional to the complex conjugate of the incident field. Therefore there exists a 
phase conjugate reflection. 

2. Scattering 

Phase conjugation by stimulated scattering was first reported by Zel'dovich, 
Popovichev, Ragul'ski and Faizullov in 1972 [Ref. 1]. As stated by Hellwarth 
[Ref. 5: p. 171]: 

In essence, the nonlinear polarization density, mediated by some driven (Raman-, 
Brillouin-, etc.) active excitation of the optical medium, couples normal-mode 
solutions of the linear Maxwell equation so as to create a new set of normal 
modes in backscattering with, in effect, complex propagation constants 
representing index change and gain. 

For example Stimulated Brillouin Scattering, or SBS, is a direct result of the 
electrostrictive effect seen in some materials. This effect is seen when the polarizability 
of a material is a function of pressure. In these materials it is possible for light to be 
scattered by a sound wave allowing the coupling of a pair of light waves to a pressure 
wave. Another type of scattering sometimes responsible for phase conjugation is 
Stimulated Raman Scattering or SRS. SRS is the scattering of light from molecular 
vibrations or rotations in a gas or liquid, or longitudinal-optical phonons in a solid. In 
SRS light waves are coupled because of polarizability changes with a molecular 
coordinate. Of particular importance to this research is Stimulated Photorefractive 







Scattering in barium titanate. In SPS it is the photorefractive effect which is 
responsible for the production of a periodic index of refraction grating; this grating 
then Bragg scatters light in the backward direction which comprises the phase 
conjugate beam. 

The mathematics involved in all three types of scattering is similar and can be 
treated together . To this end I outline the work of Hellwarth [Ref. 5: pp. 177-180] and 
discuss the existence of the phase conjugate wave assuming scalar or symmetric 
scattering with a monochromatic, multimode, incident wave in a wave guide. The 
results are largely independent of the nature of the wave guide provided that a specific 
TEM mode that propagates in the forward direction also propagates in the backward 
direction. 

A wave Ej(r) may be expressed as: 

Ejfr) * £ A m e m exp[i(k m z - pz)/2], (2.18) 

where p is the attenuation coefficient, £ m is the normalized transverse-mode pattern, 
and the propagation constants may vary with wave angular frequency, that is k m - 

U v >* 

The electric field of the backscattered wave may then be written as 
[Ref. 5: p.180] 

E 2 " \ B n e n exp[-ik n z + pz/2 ” Yz/2] , (2.19) 

where k n * k n (tD), and y is used to balance any nonlinear terms from the fields 
interaction with the medium. If y*0 in the above equation, then this is the solution 
to the linear Maxwell equations for a wave in a wave guide. 

If in Equation 2.19 y*o. then there is one configuration that allows the y of 
one of these modes to have a positive reai part significantly greater than that of any 
other. This means that a single wave (reflection) may dominate the backscattered 
waves if the necessary conditions are met. The desired result is that 



where n is a constant. Substituting Equations 2.18 and 2.19 into Equation 2.20 gives 
the polarization density necessary to balance the extra terms in Maxwell's equations 
that come from the last term in Equation 2.19 . According to Hellwarth [Ref. 9] the 
necessary conditions are satisfied if 




1-pr E H dxd y c V*V'n K mijn A i A l B » 


( 2 . 21 ) 


where 


mijn 


j^(l/L) exp[AkLz]dz, 


( 2 . 22 ) 


and 

Ak - k rn (( 0 ) + k.(v).k j (v)-k n (<D). (2.23) 

In Equation 2.21 I have assumed that the scattering considered is the pure scalar or 
symmetric type, reducing the fourth rank susceptibility tensor to a scalar quantity. 
This quantity, and all other constants, I have absorbed into t). 

The necessary conditions may be determined for special cases, either by exact 
solution or perturbation methods developed by Hellwarth. These solutions indicate 
that a backscattered solution that is nearly phase conjugate predominates because it 
has gain nearly twice that of any other mode, provided that [Ref. 5: p. 177]: 

a. the interaction length is not too long, 

b. the Stokes shift (v - to) is not too large, 

c. the number of guide modes falls within a certain large range, and 

d. the total number of guided modes is not too large. 

The solutions that led to these conditions, presented by Hellwarth in 1978 [Ref. 9], 
unfortunately Jo not rake into account competing nonlinear e:Tects. and therefore ire 
good only to the extent that they predict the major effects observed in experimentation. 
It is observed that the backscattered wave is not always an exact phase conjugate of 
the input wave, possibly due to some of these competing effects. These analyses apply 
equally as well to SPS, SBS and SRS [Ref 3]. 





C. DISTORTION CORRECTION BY PHASE CONJUGATION 

Although not immediately applicable to this investigation, the ability of a phase 

conjugator to correct for phase distortions is one of the primary aspects driving 

research in this area. As a practical point, this ability is an easy way to determine that 

any given observed beam is a phase conjugate of another. 

It was shown above that, at least for the two cases considered, a wave can be 

generated propagating directionally opposite to an input wave, with it's envelope 

function proportional to the complex conjugate of the input beam's envelope function. 

This phase conjugate wave may be considered a separate wave that must satisfy the 

same wave equation as the input wave. This means that if the wavefront of the input 

wave is distorted in some way, then the phase conjugate wave will propagate 

backwards along the same path and its envelope function will remain everywhere the 

complex conjugate of the input wave's; that is, its wave fronts will coincide with the 

input wave fronts at every point. Thus a distorted wave front will be restored to its 

original shape as the phase conjugate wave traverses the distorting medium in the 

• 

reverse direction. The proof is quite simple, as shown by Fisher and Yariv [Ref. 5: p. 
17], and is accomplished by complex conjugating the wave equation for the forward 
going wave. The result is the equation for the conjugate wave propagating in the 
opposite direction. So a wave traveling in the opposite direction of the incident wave, 
with an envelope function the complex conjugate of such a wave, satisfies the same 
wave equation as the incident wave. This means that the backward going wave will 
remain everywhere the complex conjugate of the incident wave, and therefore a 
distorted plane wave will be restored to its original configuration when it is passed back 
along the same distorting path. This distortion correction capability is the most visible 
and fascinating effect attributable to phase conjugation at the present time. 

An incomplete list of some applications of phase conjugation with emphasis on 
those currently using barium titanate as the phase-conjugator is [Ref. 5,10,11]: 

a. Holography 

b. Distortion correction 

c. Laser oscillators 

d. Interferometry 

e. Laser beamsteering 

f. Associative memory investigations 

g. Image amplification 

h. Laser coupling 








i. Image addition, subtraction and enhancement 

j. Optical bistability 

k. Automatic pointing and tracking. 

D. PHOTOREFRACTORS AND FOUR WAVE MIXING 

In photorefractive materials the index of refraction changes with incident optical 
energy. This effect was first discovered in the 1960's [Ref. 12,13] and has since led to 
the development of optical recording sensitivities comparable to that of the silver halide 
emulsions used in photography [Ref. 14]. The process involved in the effect is 
summarized by Feinberg [Ref. 5: p. 418] as follows: 

a. Light causes charge to migrate and separate in a crystalline material. 

b. The separation of charge produces a strong electrostatic field. 

c. The electrostatic field causes a change in the refractive index of the crystal by 
the linear electro-optic effect. 2 

In the photorefractive effect the important factor in the change in index of 
refraction is not the intensity of the incident light, but the total energy of.the incident 
radiation; the power only determines the speed of the reaction. Also because of the 
type of effect concerned, the actual result is dependent on the relative intensity of the 
incident radiation. 

Consider a thin beam of light incident on a photorefractive crystal that is 
otherwise dark. The light will liberate charge carriers from the lighted area of the 
crystal and these carriers will then settle back into a place within the crystal structure. 
If the carriers settle back in the area of the crystal that is lighted they will again be 
liberated. If, however, the carriers land in a dark area of the crystal, they will stay 
where they land due to the small dark conductivity of the crystal. Eventually there will 
be few charge carriers in the lighted area and many in the dark area. This 
disproportionate distribution creates an electric field between the dark area and the 
lighted area that creates a change in the index of refraction of the crystal through the 
linear electro-optic effect. 

This process of migration has been aptly described by R. A. Fisher as the 
flyswatter effect [Ref. 16j. It is as if the charge carriers were ilies at a picnic table and 
the light was a flyswatter suspended from above, slowly swinging across the table top. 

2 The linear electro-optic effect is the change in the index of refraction caused by 
(and proportional to) an applied dc electric field. This effect is seen in crystals that do 
not exhibit inversion symmetry, such as barium titanate in the tetragonal phase 
[Ref. 15: p. 275]. 





Each time a fly lands on the table within the area of the flyswatter's reach, he will 
immediately fly up as the swatter returns toward him. If, however, the fly lands on the 
table in an area away from the swatter he will be left alone and will stay on the potato 
salad. In time therefore, there will be no flies on the table in the area where the 
flyswatter can reach and many where it cannot. This process effectively pumps the 
flies to regions on the table which are away from flyswatter. 

The model of charge transport described above is known as the hopping model 
and was first proposed by Chen to explain the change in refractive indices in LiNb0 3 
and LiTa0 3 [Ref. 17]. There also exists a diffusion model proposed by Amodei 
[Ref. 18] that may dominate when the spacing between light and dark areas within the 
crystal is less than one micrometer. 

From the flyswatter analogy it is clear that what is important is the optical 
intensity relative to the ambient. This entire explanation assumes the existence of 
charge carriers in a crystalline material. In most materials exhibiting the 
photorefractive effect the origin of the charges is unknown, but it is assumed that they 
inhabit low-lying traps formed by impurity sites in the crystal. For example, in 
BaTi0 3 , Fe 2+ and Fe 3 " 1 ' impurities may create the charge carriers as shown by Klein 
and Schwartz [Ref. 19]. It has also been proposed, and there is extensive evidence to 
support the theory, that oxygen vacancies are the donors and barium vacancies are the 
acceptors [Ref. 20]. Whatever the dominant species is, it is assumed that there are a 
large number of empty receptor sites in the material available to the liberated charge 
carriers. 

If, instead of a single beam, two beams of identical wavelength and equal 
intensities (one called a pump wave and one called a probe wave) intersect within a 
photorefractive crystal, the interference pattern they form will cause a periodic optical 
field to be set up in the region of intersection, resulting in a periodic charge separation. 
Through the electro-optic effect, a periodic variation in index of refraction will result. 
This index grating is then able to Bragg scatter a pump wave in the opposite direction 
to produce a phase conjugate beam. The formation of the index grating is shown in 
F'gure 2.2, adapted from Feinberg [Ref. 5: p. 422]. 

It is important to note that photorefractive phenomena are not merely four-wave 
mixing phenomena. Four-wave mixing is a consequence of an index change at some 
given point at some given time (x,t). Photorefractive phenomena however, are 
dependent on the history of the medium and are not spatially local effects. That is, the 
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Figure 2.2 Index grating formation due to the photorefractive effect [Ref. 5). 

effects at (x,t) are a consequence of actions at (x-6x, t-5t). This is clearly seen in 
Figure 2.2. 

The top curve of Figure 2.2 shows the intensity as a function of the position 
along the path of the beam. The second curve is the charge density due to this 
intensity variation. The electric field intensity is presented in the third curve (note the 
important V 4 period shift relative to the charge). Finally the variation of the actual 
index of refraction within the crystal is shown in the last curve. 

The spatial shift between the change in index of refraction and the variation in 
intensity due to the interference pattern (one fourth of the grating period) is an 
important effect. This happens because the peak of the charge distribution is a center 
or' right-left symmetry, and therefore mere can be no electric Held at that point. This 
shift leads to the ability to transfer power from one beam to another, a process known 
as two-beam coupling. Two-beam coupling, first reported by Staebier and Amodei 
[Ref. 21], is due to the scattering of one beam into another by Bragg diffraction as they 
pass through the grating they have mutually produced. The result is that one beam 
experiences gain at the expense of the other. 












One interesting effect of four-wave mixing in photorefractors is seen in the phase 
conjugate mirror with continuous wave gain. This was first demonstrated by Feinberg 
in 1980 [Ref. 22] using the configuration shown in Figure 2.3 . 
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Figure 2.3 Formation of an optical cavity with a phase conjugate mirror [Ref. 5j. 

In this configuration a phase conjugate mirror (PCM) and an ordinary mirror (Ml) 
can form an optical cavity in which a beam will oscillate. One of the more dramatic 
demonstrations of this type of oscillation is reported by Feinberg [Ref. 5: p. 436| in a 
configuration using an ordinary kitchen spatula as one end of the cavity. Just as in 
any conventional laser oscillator, the oscillations can begin with noise (as little as a 
single photon) that makes one successful round trip within the cavity. With each 
additional round trip the system sees gain due to two-beam coupling and in a short 
period of time a visible CW beam may be observed between the mirror and the phasc- 
conjugator. 

E. STIMULATED PHOTOREFRACTIVE SCATTERING 

The second theory which accounts for phase conjugation in BaTiO, was 
originally proposed by Lam in 1985 [Ref. jj. Lam proposed That tne origin )f :iie 
phase conjugate beam is two wave mixing or scattering. This process of SPS begins 
with random inhomogeneities in the photorefractive medium which give rise to noise 
photons from the Rayleigh scattering of the input beam. In the words of Lam [Ref. 3j: 

The noise photons propagating along with [sic] axis of the input beam undergo 
parametric scattering, resulting in an unidirectional gain. The phase conjugate 








wave arises from the component of the scattered wave which experiences the 
maximum gain. 

Later the same year Chang and Hellwarth conclusively demonstrated SPS in BaTi0 3 
(Ref. 23]. To do this they had to index match the surfaces of the crystal to make 
accessable angles which would not produce the characteristic comer reflection (i.e. the 
CAT corner) producing DFWM. 

Later, Valley [Ref. 24] defined SPS as the selective amplification through the 
photorefractive effect of optical radiation scattered by medium inhomogeneities. Thus, 
SPS is a two-wave mixing process in which the incident beam interferes with scattered 
light to create the necessary grating through the photorefractive effect. As soon as the 
gain for an individual noise beam exceeds its loss, then that beam will become 
dominant, increasing the grating strength and hence the backscattered beam gain. This 
process will continue until a steady state value is reached. 

The mathematical equations, as derived by Lam, are of the identical form to 
those encountered in stimulated scattering in a waveguide (presented above), so the 
analysis of phase conjugation by SPS has already been addressed. 

F. BARIUM TITANATE 

1. The Crystal 

The experiments reported here were performed with one of the more popular 
photorefractive crystals capable of creating a phase conjugate reflection. Above the 
Curie temperature of 120°C, BaTi0 3 is a cubic crystal and therefore centrosymmetrical 
and non-photorefractive. Between 5°C and 120°C however, the crystal is in the polar 
phase and has a tetragonal symmetry. The axis of fourfold rotation is the optical axis 
or c-axis. In the tetragonal phase barium titanate is ferroelectric, photorefractive, 
photoconducting, birefringent, electro-optic, pyroelectric, piezoelectric, photoacoustic 
and (slightly) electrostrictive. The many possible effects make the crystal itself very 
difficult to understand, since often the results of the various effects cannot be 
separated. All resuits reported here involve BaTiO, in the tetragonal phase. 

Since 3a FiO, ;s chemically ind mechanically stabie at room temperature, it is 
used extensively in many capacities and is therefore well studied. Jona and Shirane 
[Ref. 25: pp. 108-215] report extensively on the structure, properties and the effects of 
temperature, stress and electric fields on barium titanate in the tetragonal phase. 
Although at the time of printing the photorefractive capabilities were unknown, their 
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data represent a comprehensive picture of the crystal and forms the basis for analysis 
of the solid state aspects of barium titanate. More recent but less comprehensive data, 
emphasizing the optical aspects of BaTi0 3 , have been compiled by Wemple, et. al. 
[Ref. 26], 

When used as a photorefractor, it is necessary that a single crystal of BaTi0 3 
be used and that the crystal be poled into a single domain. Poling is accomplished by 
cooling the crystal through the Curie point while an electric field on the order of 10 3 
V/cm is applied across the optical axis. This poling of the crystal results in an electro- 
optic tensor which has a very large [42J element. Since the photorefractive effect uses 
the linear electro-optic effect, it is important to maximize this [42] element. This 
maximization is accomplished by using extraordinarily polarized light incident at an 
angle between zero and 90° to the c-axis. 

2. Beam Fanning 

When a beam of extraordinarily polarized light is incident on a poled, single 
crystal of barium titanate, it is seen to asymmetrically defocus in the direction of the c- 
axis prior to the onset of phase conjugation. This process, called beam fanning, was 
first discovered by Feinberg in 1981 [Ref. 27]. Beam fanning is clearly visible in Figure 
2.4 (the optic-axis is designated by the arrow). This fanning can be explained by either 
an asymmetric change in index of refraction within the incident beam path, or by 
stimulated scattering. 

a. Deflection by Photorefraction 

Beam fanning was proposed to be due to an asymmetric change in the 
index of refraction caused by the photorefractive effect. This process is represented in 
the curves of Figure 2.5 (from Feinberg [Ref. 27] ). The top curve show's the intensity 
curve of a Gaussian beam. The second curve shows the net charge due to 
displacement of charge carriers by the incident light. The third curve show's the 
induced electric field, and the last curve shows the asymmetric change in the index of 
refraction. It is assumed that at the center, and most intense portion of the beam, the 
gradient of'he index of -efraction is greatest. Therefore *he center portion of the beam 
creates tne observed elfect. and tne fanning is induced beginning at that point. In 
other words, beam fanning begins in the center of the beam. This is one of the major 
discriminating factors and will be referred to later. 3 The second major point to consider 


3 Note that even if the intensity profile of the incident beam is not smooth, the 
top curve of Figure 2.5 represents the envelope that the beam intensity must fall 
within. 
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Figure 2.4 BaTi0 3 exhibiting SPPC. 
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about this theory is that as the diameter of the beam decreases the electric field 
increases, since the linear distance between charge distributions decreases. Since the 
photorefractive effect is proportional to the electric field, then beam fanning must 
increase with a decrease in beam diameter. This is clearly seen by examining the 
equations describing the electric field and the subsequent refractive index change. 

The electrostatic field induced by an incident beam in a photorefractor is 
given by [Ref. 27]: 




-2/c R T(xx + vy) , , , 

E(x) = —8-- exp((x 2 + y 2 ),'co 2 ], 

air/ ** u 


(2.24) 


where i is the charge of "he charge carrier, is Boitzmann s constant, I :» me 

temperature of the crystal, and <o Q is the incident beam diameter. The change in index 
of refraction in a crystal of BaTiO^ for an extraordinary ray due to this field is 
approximately given by 
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Figure 2.5 The process of beam fanning due to photorefraction 
for a single incident beam [Ref. 27| 
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where n(0) * (n g 2 cos0 + n Q 2 sin0f 1 ‘ , 0 is the angle the incident beam makes with 
the non c-a\is of the crystal, n o is the index of refraction for an ordinary ray, n e is the 
index of refraction of an extraordinary ray, and all terms that do not include the r 4 -, 
element of the electro-optic tensor have been dropped, since its value exceeds all other 
elements by an order of magnitude. 

Inspection of these equations clearly shows that the change in index of 
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In I9S5 Lam proposed that beam fanning is due to an effective nonlinear 
index of refraction which causes the crystal to behave like a diverging lens [Ref. ][ 
This nonlinear index of refraction t ~ is due to an intensity dependent llequcncy 
shift of the incident beam due to the bulk photovoltaic effect. Recently Valley 












(Ref. 24] has suggested that beam fanning is actually near forward scattering due to 
SPS, in an extension of Lam's theory. 

As in other types of scattering, SPS is the selective amplification of 
radiation which is scattered by inhomogeneities in the medium. The amplification is 
accomplished through the photorefractive effect. This is schematically diagrammed in 
Figure 2.6. 
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Figure 2.6 The process of SPS [Ref. 24]. 

In a crystal where one charge carrier dominates, the gain-length product 
per absorbed energy per volume at turn on is given by [Ref. 24] 


* r efT-entE d l ef r 
2n f [€C 0 (l+E d /E m )J ’ 


(2.26) 


vr.ere 'he eiieetr-c Pocxeis coefficient for extraordinarily polarized light [Ref. 2S|, 

k = 2k X, e is the carrier charge, T is the space charge field and refractive index grating 
decay rate, £ is the dielectric constant (which is anisotropic), and 
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where /eg is Boltzmann's constant, p. is the carrier mobility (which is also anisotropic) 
and T is the temperature of the crystal. Finally, l e ^- is the effective interaction length 
which Valiev models with 


l e ff " L 0 * exp[-<o o ,'(L sinB)]}, 


(2.28) 


where L is the length of the crystal in the direction bisecting the scattered and incident 


beam, <o o is the incident beam diameter, and 


8 - ( 9 , • 6 )/2 . 


(2.29) 


The important item to note is that the gain-length product is proportional 
to the effective interaction length given by Equation 2.28 . Therefore the amount of 
scattering due to SPS (both forward and backward) is exponentially dependent on the 
incident beam diameter in such a way that an increase in beam diameter will increase 
the scattering, and a decrease in beam diameter will decrease the scattering; this is 
opposite from that expected by the theory discussed in Section F.2.a above. Also, 
since SPS is dependent on the interaction of random noise photons with the incident 
beam, one would expect that beam fanning (near forward-scattering) could commence 
at the periphery of the beam in addition to the center. This scattering of energy should 
be more obvious at the edges of the incident beam since the intensity is lower there. 
Also, any light scattered toward the c-axis from inside the beam would take power 
from the periphery through two beam coupling, thus adding to the loss of intensity 
toward the beam edge. Considering these points, if the intensity profile of the incident 
beam were to be monitored during the fanning process, one should expect to see the 
intensity diminish from the periphery toward the center. Again this is opposite the 
effect expected from the prior theory. 

3. Self Pumped Phase Conjugation 

3eam fanning is intricately involved with self-pumped phase conjugation 
tSPPC). SPPC in BaTiOj was first discovered by Feinberg soon after the discovery of 
beam fanning [Ref. 2]. According to Feinberg, the process of beam fanning (due to an 
asymmetric change in index of refraction within the beam path) may result in a portion 
of the incident beam being bent in such a way that it is internally reflected off the 
comer of the crystal (see Figure 2.4). This reflected beam then intersects the incident 








beam causing a steady state interference pattern. The light and dark areas of the 
interference pattern cause an index grating (through the photorefractive effect) which 
then Bragg scatters the incident beam in the process of DFWM, creating a phase 
conjugate beam. Lam alternatively theorized that the origin of the phase conjugate 
beam is two-wave mixing (SPS) as described above, and that auxiliary beams were not 
essential to the process. 

Of interest also is the time for phase conjugation to begin in barium titanate. 
The time to produce self-pumped phase conjugation can vary from a few seconds to 
tens of minutes depending on the intensity of the incident beam and the configuration 
used. Because of the typically large time for the onset of phase conjugation in BaTi0 3 , 
there is extensive work presently underway to try to decrease it. Some work on 
changing crystals after growth has been reported [Ref. 20], but the present emphasis is 
on the doping of the crystals during growth [Ref. 29]. 

The exact processes involved in beam fanning and phase conjugation are still 
unresolved in the literature. In the following chapters I will examine these theories as 
they apply to observed effects and, since in the case of beam fanning there are 
mutually exclusive indications, I will show that beam fanning must be due solely to 
near forward SPS. 





III. PSIM: A SIMULATION OF THE PHOTOREFRACTIVE EFFECT 

IN BARIUM TITANATE 


A. INTRODUCTION 

To investigate the photorefractive effect in barium titanate, I designed a 
computer simulation which plotted the internal beam pattern within the crystal. This 
simulation, entitled PSIM (Photorefractive SIMulation), only accounts for the 
nonlinear photorefractive effect, not for SPS or other competing nonlinear effects. In 
so doing, I isolated the photorefractive effect from other effects and determined the 
necessary change in index of refraction across the beam path to account for beam 
fanning which is typically observed in the laboratory. The source code for PSIM is 
found in Appendix B. 

B. THEORY 

As discussed above, the photorefractive effect may be solely responsible for the 
observance of beam fanning, which is in turn responsible for self-pumped phase 
conjugation [Ref. 2]. The change in index of refraction of a given crystal, An, may be 
given by the equation [Ref. 14] 

An(z)~-( l / 2 )n 3 r efr E(z), (3.1) 

where n is the index of refraction, r e ^ is the effective Pockels coefficient and E(z) is the 
space charge electric field. 

Values for r e ^p range generally in the picometers per volt range. Of seven crystals 
used by Klein and Schwartz [Ref. 19], r e ^ ranged from 4.2 * 10‘ 12 to 12.0 x 10' 12 
meters per volt. Although these values are very small, the space charge electric field 
may be very high, on the order of 10 5 volts per meter [Ref. 5: p. 418]. Using these 
values it becomes obvious that the change in index of refraction due to the 
photorefractive effect in barium titanate may be as high as 10* 4 . Though not an exact 
number for any given crystal, this may be accepted as close to an upper bound for An 
(denoted here as An max ) for most of the crystals used. Gunter [Ref. 11: p. 226] gives 
An^x as 2.2 * 10' 5 , limited by the possible space charge field. 



From the above simple calculations it may be assumed that for the 
photorefractive effect to be solely responsible for self-pumped phase conjugation, the 
dimensions of any given crystal used must be such that a change in index of refraction 
on the order of 10' 4 is sufficient to cause the observed beam fanning. However, as 
previously mentioned, the gradient of the electric field and hence the index gradient is 
the important quantity relating to the photorefractive effect. Therefore, the beam size 
relative to the length of the crystal along the direction of propagation and the intensity 
profile must also play an important role in SPPC. These factors, in addition to the 
incident angle, were eventually incorporated into the simulation described here. The 
resulting simulation was used to determine the minimum change in index of refraction 
necessary to account for observed beam fanning in any given configuration. 

C. THE SIMULATION 
1. Program Logic 

The simulation is written in HP basic and designed to run on a Hewlett 
Packard 300 model computer with at least two megabytes of random access memory. 
This two dimensional model takes into account only the change in index of refraction 
across the incident beam path due to the photorefractive effect. The output is a 
graphic representation of the eventual steady state beam geometry. 

The model uses a brute force approach to calculate the observed effects which, 
although time consuming and memory intensive, uses a simple algorithm. The 
program uses only Snell's law and Huygen's principle to calculate the internal beam 
pattern. Input into the program is the range of the change of index of refraction 
across the input beam. An, the crystal and beam size (in arbitrary units), the point of 
entrance of the beam into the crystal, and the angle of incidence of the incident beam 
relative to the normal. The program computes the mean index of refraction between 
two input extremes which then corresponds to the normal index of refraction (n ~ 2.4). 
The important factor in the program is the difference between the two extremes, 
An max • Logically it is this entity which drives the photorefractive effect. Note that 
the simulation is not concerned with how An tnax ;s created. 

Once the necessary parameters are entered, the internal angle of the incident 
beam is calculated and each point on the wavefront is then assigned an index of 
refraction value that it can impart to the crystal at its location. These points are then 
grouped into pieces of the wavefront (one to five separate pieces are normally used). 





The assigned index value is dependent on An max and an assigned profile which may be 
varied (more will be said about this profile later). The indices of each point of the 
wave front are then propagated in a straight line through the crystal, at which time the 
value of the index of each point in the wavefront is imparted to the corresponding 
point within the crystal. After completing this process the wavefront is propagated 
along the same path, allowing each point to move at a speed proportional to the index 
at the point in the crystal in which it finds itself. After a short time the process of 
wavefront movement is stopped and the direction of each piece of the wavefront is 
calculated using Huygen’s principle. The original indices of the points within any given 
piece are then propagated in that new direction, imparting those values along the new 
path of the piece of the wavefront and the process begins again. The result of many 
iterations of this algorithm is a graphic output of the steady state internal beam 
pattern due to the photorefractive effect. 

2. Limitations 

There are several limitations to the model that must be discussed prior to the 
introduction of results. The first, and major, limitation is that of necessity each point 
on the wavefront cannot be treated individually. The wavefront must be treated in 
sections to determine an appropriate direction for a section and this makes the model 
only an approximation to the actual continuous case. The choice of the number of 
sections is primarily dependent on the chosen index profile. A profile similar to that 
shown in Figure 3.1 may be treated as a single piece since the direction of each 
individual piece will be the same as the direction of the entire wavefront. By contrast 
compare the profile of Figure 3.2. This profile appears to lend itself to the three pieces 
shown; however, in using only three pieces the effect around the extremum are 
minimized and it may be wise to consider the wave in five pieces. 

The difference in considering three as opposed to five pieces in the wavefront 
may be seen in Figures 3.3 and 3.4. In Figure 3.3 the index profile shown in Figure 3.2 
was used in three pieces. Note the distinctness of the three discrete sections. Figure 
3.4 shows the result of the same profile divided into five sections. It is still easy to 
note the five individual sections under consideration, however the output is noticeably 
different. Ideally it would be desirable to have a very large number of pieces used in 
wavefront calculations. However, the larger the number of pieces, the larger the 
number of points in the wavefront, and consequently the larger the beam size. To keep 
the beam/crystal ratio in the proper perspective it is necessary to increase the size of 
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Figure 3.1 Linear index profile considered as one piece. 

the crystal whenever increasing the size of the beam, and the maximum crystal size is 
limited by the available memory'. The machine used for this research was limited to 
two megabytes of random access memory, which allowed for a maximum crystal size of 
500 units. To keep appropriate beam/crystal ratios, the beam size used (generally) was 
between 20 and 40 units and the number of pieces in the wavefront ranged from one to 
eight, depending on the index profile. 

A second limiting factor in the use of this simulation is seen in cases u'here 
there is beam crossover. This happens when the index profile has a larger gradient on 
the side away from the c-axis than it does on the c-axis side. An example of this type 
of profile can be seen in Figure 3.3. In this case the beam intensity profile is 
considered to be Gaussian and the more intense section of 'he beam ci.e. the portion 
nearest the center) dominates over the weaker section. This ignores any consideration 
of two beam coupling which wrnuld actually occur. 

The final major limiting factor is the inability of the model to fully account for 
a piece of the wavefront which breaks aw’ay from the main portion of the incident 
beam. A case u'here this happens is seen in Figure 3.3 . In this case the index of 
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Figure 3.2 Nonlinear profile with possible divisions. 

refraction gradient across the piece is the same as when it was part of the entire beam. 
This may be a proper assumption if the beam was continuous during fanning; however, 
as can be seen in Figure 2.4 and many others presented here, there are often portions 
of the beam which leave the major portion of the incident beam entirely. Under the 
condition of a complete disconnection it is possible that the total index gradient should 
be based on An max and not on the index gradient as it was calculated when 'the 
wavefront piece was within the incident beam; this was not programmed into the 
model. There are presently no data available indicating the effect of a detached piece 
of a beam in contrast to the continuous case. This problem is circumvented when a 
profile such as that shown in Figure 3.1 is used so that only one wavefront piece is 
:ons;dered. 

D. RESULTS 

As noted above, the output of the simulation is dependent on the index profile, 
the beam/crystal ratio, the angle uf incidence and the position on the crystal where the 
incident beam enters. The criteria for success (of possible SPPC) of a given 






Figure 3.3 PSIM output for profile shown in Figure 3.2. 

configuration was taken to be the ability of a portion of the beam to reflect off the 
crystal corner in such a way as to reintersect itself. This would allow for the 
production of auxiliary beams which are a prerequisite to SPPC. Therefore, a result 
such as that seen in Figure 3.5 would be considered capable of SPPC, and one such as 
that seen in Figure 3.6 would not. 

Since the purpose of this investigation is to model beam fanning, the most 
important parameter is the index profile. To concentrate on this parameter, a 
beam.'crystai ratio and angle of incidence that matched one of the configurations seen 
in the laboratory were chosen. The configuration chosen was an often used 
configuration with a beam/crystal ratio of 0.064 and an incident angle of 40° . The 
incident beam entered the crystal at a point 3 / 5 way across the face of the crystal. This 
corresponds to a crystal observed in the laboratory of 5.0 mm on a side with a beam of 
size 0.32 mm, incident 3 mm from the front side of the crystal. 








Figure 3.4 PSIM output for profile shown in Figure 3.2 
but broken into five separate pieces. 

Modeling was run with three index profiles. The rationale behind the use of the 
profiles will be explained in Section E of this chapter. In each case the values of 
An max which caused successful SPPC due to beam fanning are presented. 

Following the theory of Feinberg [Ref. 2], the first index profile used was that 
shown in Figure 2.5. Successful SPPC from beam fanning was seen (Figure 3.8), 
however, note the beam fanning in 'he direction apposite the c-.axis: 'his is never 
observed in practice. 

The second beam profile used is shown in Figure 3.1. A typical successful run is 
shown in Figure 3.9. Although this approximates some observations of beam fanning, 
note the lack of a continued incident beam as one often (but not always) observes. 




Figure 3.5 Beam fanning resulting in successful SPPC. 

The final profile used is shown in Figure 3.7. SPPC due to beam fanning was 
exhibited as is seen in Figure 3.10. It was this configuration which exhibited results 
closest to those observed in the laboratory. 

E. DISCUSSION 

The first index profile is taken directly from Feinberg [Ref. 2]. To arrive at a 
profile such as this it is necessary to assume that charge migration in the 
photorefractcr is symmetric, resulting m a net charge distribution such as shown in 
Figure 2.5. The beam pattern predicted by the simulation is shown in Figure 3.S . The 
major problem with this profile is the fanning of the beam in the direction opposite the 
c-axis due to the tail of the index profile. In his original work Feinberg only considers 
the central region since it is the most intense part of the beam, ignoring the tails of the 
profile. As can be seen from the model, the tails of the profile cannot be totally 




Figure 3.6 Beam fanning not resulting in successful SPPC. 


ignored. To arrive at the observed effect with the index profile of Figure 2.5 it is 
necessary to assume that barium titanate exhibits asymmetric self-focusing (i.e. self- 
focusing only in the direction of the c-axis). If asymmetric self-focusing was present, 
then the (less intense) light subject to the profile of the tails would be unable to fan 
against the c-axis. To date there has been no reported evidence of self-focusing of any 


kind in BaTiO^; however, this does not preclude the possibility since it would only 
affect a weak portion of the beam and may not be obvious to the casuai uwener. 
Note that the fanning against the axis or. the left-hand side of 'he beam .v.raid be 

abated by two beam coupling and might not be readily discerned. 

An alternative explanation is that the crystal response to the beam is not 
characterized by the curves shown in Figure 2.5. Instead, all available charge may be 
totally excluded from all portions of the beam, resulting in the profiles shown in Figure 

3.11. 
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Figure 3.7 Index profile for simulation of beam fanning in BaTi0 3 . 

In an efFort to approximate an index profile which would exhibit the necessary 
behavior, the profile shown in Figure 3.1 was used. This beam profile was arrived at 
by making the important assumption that there is a total exclusion of charge carriers 
within the illuminated section of the crystal. This is an easy profile to work with since 
the outcome is independent of the number of pieces considered in the wavefront, as 
previously mentioned. The result of the simulation with this profile is shown in Figure 
3.9. However, there are two major problems with this profile. The first one is the lack 
of extension of the incident beam as mentioned above. Note that in most of the 
photographs presented in Appendix A there is an extension of the incident beam. This 
is almost universally, though not exclusively, true (see Figure A.5 for example). 
Second, it is unrealistic in its sharp transitions, dot;; u iiese preniems .re - incu m - :i 
the introduction of the third index profile. 

The final index profile (Figure 3.7) resulted in a beam geometry which exhibited 
all of the requisite properties. This profile was arrived at by the rounding of the 
discontinuities of the prior profile; specifically, a sine function was used to arrive at the 
configuration shown in Figure 3.7. The results of the simulation are presented in 







F igure 3 S PSIM output using the index profile of Figure 2.5 

Figure 3 10 To arrr.e at this profile it '.vould be necessary for diffusion of the charge 
carriers almost completely out oi the incident beam path. The eventual diffusion of 
charge carriers completely out of the path may actually be expected since any light of 
intensity above ambient will cause any charge curriers to eventually be totally 
liberated 4 The graphs of F igure 3 11 show the necessary charge, electric field and 
intensity profile ot the incident beam necessary to arrive at an index profile similar to 


appearance of the continued beam. 


4 I em'serg shows however, that the electric field may be proportional to tic 
gradient ct the incident beam inter,'.;''. Ref 2~i 
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To arrive at the configuration of Figure 3.9 a value of An max m 0.03 was used. 
In Figure 3.10 - 0.02. Regardless of the profile used, these values of An max 

are characteristic of the value necessary to account for the observed beam fanning in 
BaTiOj. For successful SPPC with the configuration chosen, An max must be greater 
than 10' 3 . Reducing the beam size will obviously reduce the necessary value of 
Arimax- b ut values On the order of 10‘ 2 to 10‘ 3 are typical. 

Thus, the simulation has predicted the necessary change in index of refraction 
across the beam path, and the necessary profile of that change. In Chaper V the 
results of the necessary value for An max will be compared to those possible for the 
crystal. Also, the physical implications of the necessary profile of the index of 
refraction within the beam path, as determined by the simulation, will be discussed. 
The comparison of the information derived from the simulation to that experimentally 
derived from the crystal, will show that the necessary conditions cannot be met in 
BaTiO^. 




























IV. EXPERIMENTS 


A. INTRODUCTION 

In order to investigate beam fanning in BaTiOj, an experiment was designed to 
measure the influence of a change in the incident beam diameter on the fanning 
process. If fanning increases with a decrease in beam diameter, it could be concluded 
that beam fanning is a direct result of the change in index of refraction within the 
beam path, as predicted in Equation 2.25. If, however, the fanning decreases with a 
decrease in incident beam diameter, it could be concluded that beam fanning is actually 
forward-scattering due (probably) to SPS, since, as shown in Equation 2.28, the gain 
per unit length will decrease with a decrease in spot size. 

B. EXPERIMENTAL APPARATUS 

The apparatus used for this work is shown in Figure 4.1. The laser was a Spectra 
Physics Ar + variable power (20 mW maximum), multimode laser, model #162A-07. 
During the experimentation the laser was operated at all of the visible lines and was 
not optically isolated from the experiment. Upon leaving the laser, the beam was 
directed into a polarization rotator in order to rotate the beam from the ordinary to 
the extraordinary polarization in the crystal. To ensure that only extraordinarily 
polarized light was used, the beam was then directed through a Glan-Thompson 
polarizer. Following the polarizer was a beam splitter, to enable the observation of 
any phase-conjugate reflection, and a variable aperture of the type common in 
photographic lenses. This aperture acted as a spatial filter, reducing the amount of 
scattered light reaching the crystal. The beam was then focused onto the crystal by 
uncoaied lenses of either 34 cm or 16 cm focal length. 

The crystal was one of four single crystals of BaTi0 3 manufactured by Sanders 
Associates and poled into a single domain. Two of the crystals ion loan from the Los 
Alamos National Laboratories) measured approximately 5x5x5 mm and were 
transparent with a slight yellow tint. Two of the crystals (belonging to the Naval 
Postgraduate School) measured approximately 5 x 5 x 2.5 mm and had no noticable 
tint to them. The crystal was mounted on a rotatable platform capable of being 
translated 12 cm along the axis of the optical bench on which all of the components 
were mounted. Following the crystal, and connected to the mounting table, was a 



knife edge. The knife edge was coated with black tape to ensure that a resonating 
cavity was not formed with the crystal. Behind and to the side of the knife edge, also 
connected to the translation table, was a silicon photodiode connected to an 
electrometer measuring the short circuit current. 

C. EXPERIMENTAL PROCEDURE 

To measure the amount of beam fanning quantitatively, the knife edge was 
aligned in such a way that the beam transmitted through the crystal was blocked from 
the detector. To ensure that the blockage of the transmitted beam was sufficient, the 
table was translated rapidly along its entire range of motion and the detector output 
was monitored to ensure that there was little or no change in the recorded intensity. 
The vibration inherent in the movement and the changing beam diameter of the 
incident beam ensured the lack of significant beam fanning, so that a lack of change in 
the recorded intensity indicated that the transmitted beam was completely blocked 
from the detector for all observed beam diameters. 

Once the crystal, knife edge and detector were in position, the beam was blocked 
and the crystal was hooded with light from a 75 watt incandescent lamp located 
approximately 12 cm away. After approximately one minute the light was turned off 
and the laser beam was unblocked. The intensity of the light reaching the detector was 
then monitored and the peak and steady-state intensities were recorded. Once the 
intensity had reached a steady-state value the beam was again blocked, the crystal was 
translated and flooded with white light, and the procedure was repeated. 

The flooding of the crystal with white light effectively homogenized the charge 
distribution within the crystal and ensured that one reading was not affected by the 
previous one. Special care was taken to ensure that the beam was parallel to the 
crystal table and that all components were mounted on the optical axis of the optical 
bench. The precision translation mechanism on which the crystal table was mounted 
ensured that the beam was incident on the crystal with the same center for all 
measurements. These precautions effectively compensated for any surface 
inconsistencies on '.he crystal face except at very -;maii laser beam diameters. 

Both peak and steady-state values of the fanned beam intensity were recorded; 
however, only the peak values were used for final analysis. This decision was made 
because, although the steady state values generally correlated with the peak values, the 
steady state values could be affected by numerous variables. One obvious problem 
with steady-state readings is the possibility of phase conjugation, thus reducing the 






amount of fanning recorded in the steady state. Also, due to the orientation of the 
detector, it was possible (even likely) that the beam would fan across the detector, 
eventually coming to rest past the area it can record; sometimes the fanning was 
observed to exit the face perpendicular to the face the detector was monitoring. To 
reduce all of these effects, the incident beam was originally allowed to enter at only 
small angles to the c-axis (~2° - 5°) and, although alignment was tedious, consistent 
results were obtained for various crystals, lens focal lengths, laser lines and incident 
angles. Eventually larger incident angles were used (~ 20° - 40°) with similarly 
consistent results. 

To ensure that the observed effect is not influenced by scattering initiated by 
small (e.g. molecular size) crystal inhomogeneities, one configuration was chosen in 
which there was moderate fanning and the fanned intensity was recorded for each of 
the available visible lines of the laser. The result was a mean ratio of incident intensity 
to fanned intensity that was constant across all available wavelengths, thus indicating 
no resonances and no dependence of scattering on wavelength at the available laser 
lines. 

After the data were collected, the beam diameter at each point was measured 
with an Aeronca Electronics Laser Blade and the direct laser intensity was recorded. 
Here, the beam diameter is defined as twice the distance from the center of the beam to 
the point of 10% intensity. For analysis purposes the beam diameter could be divided 
by the cosine of the angle of incidence so that the beam diameter would be that 
actually along the y-axis of the crystal. 

Measurements were also made by placing the detector above the crystal as the 
incident beam diameter was changed. This gave an indication of the total amount of 
scattering, not just asymmetric scattering. In this case, the front surface of the crystal 
was shielded from the detector with black tape to keep specular reflection off of surface 
irregularities from influencing the measurements. Transmission measurements were 
also made by replacing the knife edge with an aperture, thus allowing only the 
transmitted portion of the incident beam to reach the detector. 

D. RESULTS 

Typical results of the fanning measurements are shown in Figure 4.2. In this 
case the crystal was at an angle of < 5° so the beam diameter need not be adjusted for 
cos0. All intensity recordings have been normalized for comparison purposes. The 
shape of the curves is insensitive to the laser line used, beam angle, choice of crystal, 
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and the divergence or convergence of the incident beam (i.e. whether the crystal is 
placed in front of or behind the focal point of the focusing lens). 

Figure 4.3 is a plot of the intensity seen from above the crystal for various beam 
diameters for converging light. There are two important things to note here. First, 
there is less scattering for smaller beam diameters and second, the sharp discontinuity 
at about three micron beam diameter. This discontinuity is noticable as the crystal 
goes through the focal point of the focusing lens. Figure 4.4 is a plot of the scattering 
as viewed from above the crystal for diverging light. Note the initial decrease in the 
scattering within the first four microns. This is only apparent with diverging light and 
is seen in all crystals I used. Note also that there is a much broader dynamic range in 
the intensity for diverging light than for converging light. 

Although all of the scattering data are generally insensitive to which crystal is 
used (i.e. the general shape of all the curves for a given situation is the same for all 
crystals), the same cannot be said for the transmitted beam. Figures 4.5 and 4.6 show 
the transmitted intensity for all four crystals plotted together. Crystal #2 is the only 
one which does not exhibit SPPC and is the only crystal in which the transmitted beam 
intensity increases with beam diameter as a general trend. However, note that all of 
the curves exhibit some inflection point around six micron diameter for diverging light. 
This inflection is not a consistent maximum or minimum; it depends on the crystal. 
Presently I cannot explain the cause of these dissimilar curves. Generally it may be 
noted that the transmitted intensities map the same characteristics as the scattered 
intensity, but inverted, as would be expected. All of these various features will be 
addressed later, but the important aspect of all of the data is obvious; increased beam 
diameter results in increased beam fanning. 

Although these data alone are a strong indication of the lack of index change 
within the beam path due to the photorefractive effect, there exists even stronger 
evidence. By removing the knife edge and detector from behind the crystal, the image 
of the transmitted spot could be observed on a screen approximately three meters 
v.vav. Figure <hows -he image of the spot, which was projected directly onto 

photographic .:1m a r . diiTerent times after the beam was allowed to enter the crystal, 
and prior to the onset of phase conjugation. Figure 4.7a is an image of the beam 
immediately after turn on; notice the spot is generally undistorted. Figure 4.7b is a 
photograph between turn on and steady-state, and Figure 4.7c is a photograph at 
steady-state. By viewing these images it is clear that the energy' of the beam is first 
extracted from the outside of the beam and not from the inside. 
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In addition to the above measurements, the index of refraction of the crystal was 
estimated by measuring Brewester's angle. At X * 514 nm the published value for the 
extraordinary index of refraction is n e = 2.424 [Ref. 22], however, my measurements 
consistently resulted in n e = 2.26 ±0.02. This large discrepancy is still unexplained. It 
is also interesting to note that, although there was a distinct minimum at Brewester's 
angle, the reflection did not disappear entirely. The minimum was so slight that it was 
quite difficult to find, and eventually a photodetector had to be used to accurately 
locate it. It is possible that both of these observations may be explained by the 
presence of impurities on the surface of the crystal. Specifically, an oil film from the 
hands of the experimenter may be at fault. Additionally, there is the possibility that 
heating at the surface of the crystal due to the high energy density of the incident beam 
may be responsible. Since measurements by others [Ref. 30] of the index of refraction 
of BaTi0 3 using Brewester's angle have resulted in indices consistent with the reported 
values, either one or both of these explanations may account for the observed 
anomalies. 
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Figure 4.7a Photograph of projected spot immediately upon turning on the incident beam. 



































V. ANALYSIS 


A. DISCUSSION OF RESULTS 

The simulation in Chapter III and the experiments described in Chapter IV 
clearly show that asymmetric self-defocusing cannot be due solely to a change in index 
of refraction along the beam path within the crystal. Even with the changes of the 
profile for the change in index of refraction proposed in Chapter III, the necessary 
value of An 1Tiax is too large according to the published maximum values for BaTiOj. 
Inspection of the necessary values of An max from the simulation shows that values on 
the order of 10' 2 are necessary to account for the observed fanning within the crystal; 
however, previously calculated values place an upper bound of ~ 10' 4 , and a value of 
~ 10‘ 5 is probably more reasonable 5 [Ref. 11). 

This discrepancy may be accounted for only by assuming that the index gradient 
is limited to a very small portion (in some cases as small as V. Q the size of the input 
beam) of the center of the beam, or by assuming that the observed effect is a result of 
another process. It is impossible to justify the first assumption and still achieve the 
necessary index profile, since charge carriers must be liberated from the entire region 
where the brightness of the incident beam exceeds ambient. Therefore, by Laplace's 
equation, the gradient of the electric field within the beam path must be uniform. This 
leaves only the possibility that another process is responsible for beam fanning. 

The results of the experiments of Chapter IV reinforce the conclusions made 
from the simulation portion. Figure 4.2 shows that beam fanning decreases with a 
decrease in beam diameter. From Equation 2.25 it is clear that this behavior cannot be 
explained by a change in the index of refraction within the beam path due to the 
photorefractive effect. If this type of change were responsible for beam fanning, the 
the amount of fanning would have to increase with a decrease in beam diameter. 
\visiit::r.ailv. is hour, n Figure 'earn running begins an 'he pcr.phcrv md ‘he 
.:gh: :> progressing .rom T.e >ut -,u*e v. *nc 'eair. to warn ::ie renter. Noting 

the index profile as shown in Figure 2.2, this is impossible if the fanning is due to an 
asymmetric change in index of refraction across the beam path, as was originally 
predicted. 

5 Feinberg does, however, calculate a possible maximum value on the order of 
10° [Ref. 2" T ). 











One may argue that all of the observed qualities of beam fanning may be 
explained by assuming "small-scale" self-focusing of the incident beam within the 
crystal. That is, the incident beam may experience a breakup into spots which are of 
diameters small enough to create index gradients of the magnitude shown necessary by 
the simulation in Chapter III. This type of self-focusing has been known for over two 
decades, and is common in materials whose dielectric constant increases with field 
intensity [Ref. 31]. If small-scale self-focusing occurred in BaTi0 3 , it could explain 
many of the effects that are observed during SPPC (e.g. the presence of small filaments, 
beam fanning, etc.). However, there are serious problems with a small-scale self- 
focusing theory. First, the breakup of optical beams due to self-focusing has been 
shown to be spatially periodic, with the periodicity being dependent upon intensity 
[Ref. 32]. As is demonstrated in Figure 4.7, there is no indication of spatial periodicity 
in the beam as it is transmitted through a BaTi0 3 crystal. Also, even in good self- 
focusing mediums at very high beam intensities (~ 10 7 W/cm 2 ), it is necessary for the 
beam to have a long path length (on the order of tens of centimeters) within the 
medium before self-focusing is observed [Ref. 33]. Even with a very large change in the 
dielectric constant, it is unlikely that self-focusing could occur within a medium of only 
5 mm in length. Finally, there is no evidence in the literature that the dielectric 
constant of BaTi0 3 changes significantly with the intensity of the incident light 
(although it does change with an intensity gradient, as discussed in Chapter II). It is 
interesting to note however, that the dielectric constant of BaTi0 3 does change with 
temperature [Ref. 25], and this may cause self-focusing at very high powers and very 
small beam diameters, given a long enough path within the medium. 

From the evidence presented here it follows that asymmetric self-defocusing in 
BaTi0 3 cannot be due to an asymmetric change in index of refraction along the path 
of the incident beam due to the photorefractive effect. Attributing this phenomenon to 
some other process however, is not so easy. 


B. SCATTERING 


'eu .r. Chanter II. 'tie ;mv other thecrv that s current:*-' ivaiiab.e '0 


explain asymmetric self-defocusing in BaTi0 3 , besides the theory just discussed, is the 
theory involving scattering (SPS) [Ref. 24]. Although I have not conclusively 
demonstrated that beam fanning in BaTi0 3 is due to SPS, there are a number of 
observations which indicate that this type of scattering maybe the origin of the effect. 
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As may be seen in Equation 2.28, if near-forward scattering is the origin of beam 
fanning, then there should be a decrease in fanning with a decrease in incident beam 
diameter, an observation borne out in our experiments. Also, the diversion of energy 
from the incident beam could progress from the periphery toward the center, since 
noise photons originating from outside the beam must interact with the first portion of 
the beam they encounter (Figure 4.7). 

Another interesting and important discovery is that although the shape of the 
fanned intensity curves is similar regardless of convergence or divergence of the 
incident beam (in front of or behind the focal point of the focusing lens), the absolute 
intensity of the scattered light is always observed to be greater for the diverging case. 
This is not perplexing when considering near-forward scattering. One would expect 
more scattering from diverging light than converging because the size of the beam 
increases as it progresses through the crystal for diverging light, thus increasing the 
amount of scattering. Also it is noted that individual intensity readings are not 
consistent, even if there is no no movement of the apparatus between readings. This is 
expected since the scattering process begins with noise photons, and no two situations 
can be exactly alike, even for similar configurations. This inability to exactly reproduce 
any given measurement suggests that the scattering may begin with random index 
changes within the crystal, caused by the stochastic distribution of charge carriers. 

Finally, if stimulated scattering is responsible for beam fanning, then one would 
not expect continuous fanning, but fanning only along paths from which an incident 
photon interacted with the incident beam, thus creating a set of filaments emanating 
from the incident beam. A close examination of any of the photographs of Appendix 
A shows that this is indeed the case. Fanning is a set of discrete small beams all 
emanating from the incident beam. We must now discuss curved appearance of the 
dominant-scattered beam during self-pumped phase conjugation, a phenomenon not 
explained by SPS. 

A very close examination of many photographs of SPPC (e.g. Figure A.5) reveals 
that beams which appear to curve often have a series of discontinuities within them. 
Generally, x filament. or set of filaments, leaves the incident beam along r he direction 
of the c-axis. After a certain distance (usually about halfway to the crystal edge) these 
filaments enter a region where they become fuzzy, and often two become 
indistinguishable. The filaments exit this fuzzy region at an angle different than their 
entrance angle. This type of behavior is explainable by assuming that within this fuzzy 








region there is stimulated scattering due to SPS. Forward scattered beams enter into a 
region where scattering occurs due to photons arriving from a different angle. The 
beams which continue on at a new angle are the forward scattered beams from this 
second scattering. It is therefore a series of this type of scattering discontinuities which 
gives the characteristic curved appearance of a filament during SPPC. However, this is 
not always the case. 

There exist configurations where this series of discontinuities does not appear to 
account for the curved filaments (i.e. Figure A.6). In a majority of these cases a close 
examination reveals that a large filament originates at the incident beam or a series of 
filaments leave the incident beam very close together. The curved appearance comes 
from the coupling of energy to the filament, or portion of the filament, in the direction 
of the c-axis as the filaments progress through the crystal. Therefore the beam doesn’t 
actually curve; it simply appears to do so because it is losing energy from one side and 
transferring it to the other. The process is actually a form of self-focusing due to two 
beam coupling. 

I have observed a very few instances where it appears that neither of these 
processes can fully account for the curved appearance of the fanned filaments. In 
these cases one would expect the effect to be due to true beam fanning from the 
change in index of refraction within the filament path. Although the simulation PSIM 
shows that it is impossible for this phenomenon to fully account for the observed 
fanning effect, it would be surprising not to find some evidence of a change in index of 
refraction on the order of 10' 5 across a beam of diameter less than a few microns. 

C. SELF PUMPED PHASE CONJUGATION 

From the discussion of the previous section we may conclude that a change in 
index of refraction within the incident beam path is not responsible for asymmetric self 
defocusing in BaTiOj. It is probable that SPS is responsible for the effect. Therefore 
one may say that the true origin of SPPC in BaTiO^ is scattering, probably SPS. The 
cause of the final phase conjugate reflection is not as clear. 

Phase conjugation by stimulated scattering in 3a HO, has been demonstrated 
[Ref. 23], and recently it has been shown how backscattering may dominate forward 
scattering in SPS [Ref. 24]. Additionally, the investigations reported here indicate that 
scattering is responsible for the phenomenon of beam fanning. Therefore one might 
incorrectly conclude that phase conjugation in BaTiOj is due to stimulated scattering. 
With the onset of backward scattering in the crystal there begins a competition of the 



scattered waves. As the modes begin to build, a mode showing high gain can usurp 
power from one showing less gain. In this way only the modes showing the highest 
gain eventually survive. It is this competition between configurations having high gain 
which is responsible for the amplitude oscillations often seen in a phase conjugate 
reflections (see [Ref. 34,35] or Appendix A). 

If a filament, arising from near forward scattering, has a configuration such that 
it can reflect off the comer of the crystal (i.e. a CAT comer) to create a situation 
where degenerate four wave mixing (DFWM) is possible in addition to SPS, then one 
would expect that filament to have a gain advantage. This then would be the filament 
to eventually win the competition. So it actually would be expected that both DFWM 
and SPS should contribute to SPPC. 

A review of the conditions of the only published report of phase conjugation due 
solely to SPS [Ref. 23] shows that special precautions were taken to ensure maximum 
coupling. Specifically, an index matching liquid was used which preserved as much of 
the phase conjugate return as possible. Additionally the crystal was surrounded by a 
glass cuvette (containing the index matching liquid) which increased the amount of 
scattered light available for SPS to begin. Even with these advantages, a phase 
conjugate reflection of only ~ 10% was obtained. This is quite small compared to the 
30% to 50% usually seen from SPPC in BaTiOj in air. This indicates that DFWM 
from the internal comer reflection is the most significant contributor to the phase 
conjugate return observed in SPPC, as was originally proposed [Ref. 2]. 






VI. CONCLUSIONS 


Experimental evidence does not support the theory that asymmetric self 
defocusing in BaTiC> 3 is a result of an asymmetric change in the index of refraction 
within the path of the incident beam. Instead, forward stimulated scattering, probably 
Stimulated Photorefractive Scattering (SPS), is found to be a more likely explanation. 

Simulation has shown that for an index change to be responsible for the 
phenomenon of beam fanning, the necessary change in the index of refraction across 
the beam must exceed the theoretical maximum, based on the maximum possible space 
charge field. A change in the index of refraction on the order of 10' 2 is necessary to 
account for the observed effects, whereas the accepted maximum is on the order of 10‘ 4 
to 10' 5 . 

Experimental results do not support a theory which requires a change in index of 
refraction across the incident beam as the cause of beam fanning. The magnitude of 
fanning is observed to increase with an increase in beam diameter. Also, observations 
of the manner in which the beam fans show that energy is scattered out of the beam 
path beginning with the periphery and not beginning in the center. Neither of these 
observations is consistent with the results expected if a purely photorefractive effect 
were responsible for asymmetric self defocusing. 

Close observations of the process of beam fanning in BaTiOj indicate that beam 
fanning can be explained by a stimulated scattering theory. The presence of filaments 
rather than continuous fanning, scattering beginning with the periphery of the incident 
beam, and discontinuous filaments during phase conjugation are all explicable when 
beam fanning is attributed to stimulated scattering. 

Future experiments in this area should include probing the "fuzzy" regions of the 
observed filaments from above with a low power laser to search for the presence of 
additional index gratings. Also, the actual change in the index of refraction across the 
incident oeam path within the crystal remains to be determined experimentally. 
Additionally, further research into the size, type and spacing of scattering centers 
within the crystal is necessary to fully explain the observed difference in scattering of 
diverging and converging light. 


65 




Mad 








APPENDIX A 

INTERNAL BEAM PRODUCTION 




V 

B 


1 


f§ 


1. INTRODUCTION 

The production of internal beams in BaTi0 3 has, to my knowledge, never been 
systematically investigated. A search of the available literature indicates that the cause 
for this oversight may be that, in general, millimeter sized, collimated beams have been 
used. In all available photographs of BaTiO^ used as a SPPCM, the beam has not 
been focused onto the crystal. If the beam is not focused, then the diameter of the 
beam upon entering the crystal is at best the same as when it leaves the laser. It 
appears that when a beam is not focused onto the crystal, auxiliary beams are 
eliminated or obscured by the strong input beam since its diameter is large compared 
to the dimensions of the crystal (sometimes being more than half the size of the 
crystal). However, if the input beam is focused, the beam diameter is small compared 
to the crystal and the immediate effects of the input beam stay localized within the 
crystal, thus allowing the formation and observation of auxiliary beams. There is also 
an immense increase in the energy density of the beam within the crystal upon 
focusing. 

It was using a focused beam that we first discovered the ability of barium 
titanate to support internal auxiliary beam patterns. Extensive investigation by A. V. 
Nowak and myself at Los Alamos National Laboratories indicated that there is a 
plethora of internal beams which may be generated by the nonlinear effects of barium 
titanate during phase conjugation. It appears that phase conjugation is necessary to 
sustain these beams since no auxiliary beams were observed in the absence of phase 
conjugation. In many cases it seems that the beam is initiated by the reflection of the 
phase-conjugate beam off of the face of the crystal as it exits, but this is not always 
true. 

'•Vhat follows !n this chapter is an attempt to report, categorize. and :\r!;ar. 
many of the observed internal beam patterns ana anomalous effects. It :s :mror.m: 'o 

note that often it is difficult or impossible to recreate some of the observed effects. It 
appears that not only position and beam profile are important, but also the short term 
history of the crystal. This means that possibly one effect may not be observable 
unless another has (or has not) immediately preceded it. 








2. EXPERIMENTAL APPARATUS 


All investigations presented in this chapter were performed at Los Alamos 
National Laboratory on an air floated optical bench in free air that was maintained at 
a temperature of 22 ± 1°C. The apparatus used is shown schematically in Figure A.I. 

During all experimentation the light source was a Spectra Physics 2020-03 argon 
ion laser, prism-tuned to 514.5 nm and operated in a single longitudinal mode by 
means of an intercavity etalon. The beam at the laser had a diameter of 0.88 mm at 
the 1 / e of peak intensity. 

After leaving the laser, the beam was directed through a half-wave plate that 
served as a beam attenuator when necessary. Following the half-wave plate was an 
apparatus for isolation of the laser from the phase-conjugate beam. The apparatus 
used was a faradav rotator, manufactured by Optics For Research Inc., consisting of a 
25mm long rod of Hoya FR5 glass surrounded by a set of permanent magnets. The 
isolator was set for optimal rotation of the plane of polarization of the beam at 45° 
and placed between two polarizers with axes oriented at 45° to each other. This 
configuration gave a backward attenuation of approximately 500. Without the isolator 
the phase-conjugate beam entered the laser cavity and caused large fluctuations in 
power output as well as mode hopping. The polarization of the beam upon leaving the 
isolator was such that it would enter the crystal extraordinarily polarized, which is 
required for the formation of a SPPCM. 

After passing through the isolator the beam was expanded by a factor of 6.25 
with an expanding telescope. This allowed the insertion of transparencies into the 
beam when desired, enabling the verification that an observed beam is a phase- 
conjugate beam. Also included in the telescope apparatus was a 0.05 mm pinhole that 
could be removed at will. This pinhole provided the ability to make observations with 
a Gaussian beam (pinhole inserted) or with what w r as eventually termed a mottled 
Gaussian beam (without the pinhole). The mottled Gaussian was so named because, 
although the beam without the pinhole was generally Gaussian, it had a mottled 

> •n»'i —io<- > 

:w:nc oore*, tne 'earn was iirecteu through a beamsplitter with o7 ,J . 

transmission at 514 nm. This beamsplitter not only provided the means for viewing the 
phase-corrugate signal as it returned along the incident path, but also provided a 
reference beam that was used to look for frequency shifts of the phase-conjugate beam. 
To iOck for a frequency shift the phase-conjugate return was reflected from the 
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Figure A.l Experimental Apparatus. 













beamsplitter and redirected to interfere with a portion of the incident beam reflected 
from the front side of the beamsplitter. The resulting interference pattern was then 
observed for evidence of movement. 

Upon passage through the beamsplitter, the beam, with power of 75 ±5 mW, 
was erher modified by a lens, slit, knife-edge or a combination thereof, or simply 
collimated and directed toward the crystal. As the phase-conjugate beam returned 
along the incident path, the beamsplitter reflected a portion of it and directed it to a 
large area silicon PIN photodiode reverse biased at 22 volts and terminated with a 100 
ohm resistor The output of the detector was sent to a digital millivolt meter, an 
oscilloscope and a chart recorder. When an image was placed in the beam the detector 
was removed and replaced with a white screen for viewing. W'hen frequency shift 
information was desired, the detector was replaced with a mirror and the two beams 
were directed toward a white screen for observation of the resulting interference 
pattern 

The cry stal was any one of four nearly cubic single crystal of BaTi0 3 measuring 
5±'V2 mm on a side. The sides of the crystals were polished, but it is important to 
note that the sides were not exactly the same size and therefore the sides were not 
exactly orthogonal. The crystals were mounted on a platform that could be raised, 
lowered, rotated and tilted. Observations of internal beams were made visually from 
above and were recorded by film or video tape Both the still and video cameras used a 
f * 3 <; mm Cannon FL lens mounted on a single track bellows extended to provide a 
full frame view of'he crystal 

Many of the photographs presented here were taken from the video tape directly 
off c- i television screen using a cm format still camera The desired frame of the 
video was dign/ed and placed on a hard disk for recall and photographing, the 
numbers displaced in 'he photographs are a resuit of the equipment used to do this 
final pho'ograrhs ind ire nor pertinent to the investigation When this method was 
-'ed. / a is pussmie to m-iKC negative images of the desired photographs, which 
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±3°o. In the photographs presented, except where noted otherwise, the incident beam 
enters vertically from the top and is visible only upon entering the crystal. 

3. RESULTS AND ANALYSIS OF DATA 
a. Experimental variables 

The number of variables involved in SPPC is immense, however I have tried to 
concentrate on those which have emerged as the most important. The most important 
factors in the production of vanous beams differ, but they generally may be considered 
as dependent on beam profile or the geometry of the beam-crystal interface. 

Lenses, knife-edges, slits and pinholes were used to change the profile of the 
incident beam prior to the beam entering the crystal. The effect of these modifications 
of the incident beam is especially obvious when the beam profile has a large gradient, 
such as when a knife-edge is inserted into the beam, since this creates a very clear 
delineation in the lighted and unlighted areas of the crystal and introduces new wave 
vectors into the beam. The effect of the lens to crystal distance on phase-conjugate 
reflectivity has been studied by A. V. Nowak and is still not well understood [Ref. 37], 

The geometry of the beam-crystal interface is probably the most important of 
all variables in the system. Changes in angle of less than one degree and lateral shifts 
in the beam of less than 0.5 mm can cause drastic differences in internal beams 
configurations. The following examples are presented to indicate the drastic effects 
observed with these type of changes. 

The configuration shown in Figure A.2 is a typical experimental configuration 
that produced a phase-conjugate signal from the BaTiOj crystal. The phase conjugate 
reflectivity was 33%. Note the fanned beams at the lower left comer of the crystal. 

Figure A.3 demonstrates the extreme effects made possible when changing the 
beam profile by placing a knife-edge in the incident beam. Unless otherwise indicated, 
the knife-edge always intercepted the beam on the right side and was placed behind any 
lenses Note the extreme curvature of the internal beam observed when blocking half 
the beam with a knife-edge This also reduced the phase-conjugate reflection to zero 

\ <arge effect is aiso seen when relocating the beam entrance position by a 
iarerii shift in the bej.m-cr.stal intersection As the beam is moved across the face of 
tne crystal, the phase-conjugate reflectivity will change, beginning with no phase 
ccn;uga'i;-n it one end and ending with the same resuit at the other In between the 








Figure A.2 Normal configuration for BaTiO } exhibiting SPPC. 


two extremes, the phase-conjugate reflectivity may range as high as 60% or be 
nonexistent depending on other factors. Similarly, a change in the angle of the beam 
entering the crystal will have a drastic effect on the phase-conjugate reflect!vity and 
beam pattern. It is generally unpredictable and small changes in angle ( — > ha\e 

been observed to have no effect, to destroy a pattern and eliminate all pnase-conjugate 
reflectivity, or to increase the phase-conjugate reflectivity as much as 500%. 
b Complex Auxiliary Beam Patterns 

/. The Total Internal Reflection ( TIR) Ring 

One of the most interesting effects discovered was the TIR ring Figure 
A J shows the first example, to my knowledge, of a totally internally pumped ring 
an within i BiTiO. .r <\ii. This effect '-a- iff:.eve : ^s.r.c i :e- 

.cr.>. i'f ? v.m irom :ne ;ace ot tne crystal I he beam was incident at -iff’ to the axis of 
the crystal and the knife-edge was placed on the side of the beam nearest the cfvtal 
The effect was occasionally obtained while obscuring the other half of the beam, and 
was once obtained with a vertically aliened slit A V Nowak once obtained the rme 
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Figure A.3 Effect of placing a knife edge m the incident beam. 

with just the focusing lens, indicating that the TIR ring may be dependent more on the 
availability of light scattered at a certain angle than on the incident beam profile 

In the case shown in Figure A.4 the phase-conjugate signal took about ><i 
seconds to appear The beam within the crystal slowly moved to the corner until the 
configuration of the beam arrived as shown, but without the ring. Then. slow|\. the 
ring appeared over a three second interval. As the ring appeared, the phase-conjugate 
retlectsvitv dropped from "’6' > o to 6"°. A drop in phase-conjugate refectr. uas 
observed :n all examples of ring formation for those cases where the refectr-ity us 
measured In some cases, several nearly overlapping rings were observed, as in 1 igure 
A 4 The ring in Figure A 4 was steady, fading occasionally ion the order o: rue 
er- ' • e rrunu*es• '-u* fv. r amir.a '* s i>. rot wer: 1 e ' ■■ ,■ 

: • ■ - ■ ■ , . • . .... 

iti'.g *-t t.ie ratutiun was anuer i tncr. tne ring would reappear 1 ians.at.on ac'ng 
the cr.stui on 'r.e order of «> ' mm in either direct.on did not disturb the ring. t„rhe- 
translation would destroy it f or other configurations, the ring was or..\ preset .eu ! r 
one direct.on of cr.stai motion ■ *ee Secti m J> 










Figure A.4 The TIR Ring. 

The direction of circulation within the rings wa» determined b\ extracting 
some of the light out of the ring into a second crvstal of BaTiOj After producing the 
ring m one crvstai. a drop of index matching liquid imcthvlene iodide CH,I,' dilu’ed 
1 1 with 2-propanol, n — 1 6> was placed onto the surface of a second cnstal and ’he 
second crxstal was moved into contact with the first Since the loss into the second 
.nstu. was too great to allow the continuation of the ring, the operation was recorded 
K r ’he video camera and the frame of the .nstant was of contact was stud.ed a' i .uter 
’.me St^d - . of the video revealed that the light within the ring could pr< paga'e 
counterclockwise, or both wavs simuitaneouslv depending on the part.cw.tr • n.g 
Howes er. n 'he majoritv of cases, oniv ’he counterclockwise direction was '^nc; 
f c \ ’ hear 1 '. >h ss *h.s . , .r.’er lo.kw se oug.it . 1 n .r.d - :vrrt • •• 


1 fie ring withm the cr.st j. ,s pruhuhls produced o\ tvso-Ke v . 
he’ween 'he main beam and scattered ught that finds a closed l'"p ■. ,.i 
'ctlecticns around ’he crvstul l his somewhat . ss\ MR ring takes power ’• 
..undent u cjh, which explains the drop in p’U'C-c<.n;ugate reflect. v .tv * een i 
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Figure A.6 The Diamond configuration 

sharp diamond appears in the crystal the phase-conjugate signal increased considerabh 
up to >0 times .shat it is just prior to its appearance) The sharp diamond prohibit 
begins is a specular refection if the phase-conjugate signal, winch reflects back to the 
m.cident beam b>, internal reflection Since a second point for DFWM is created, one 
•' o .id expect the observed increase in the phase-conjugate signal upon the appearance 
' f ':.c diamond 

A fascinating and unexplained sanation o; the sharp diamond is the dill use 
d.irr.ond I he Jifiuse diamond, shown m Figure A begins m a configuration simalar 
• : ",e di imm.J shown in Figure A o. but i' the crvs'ji ;s lateralis moved, a portion of 
me f;:',t '->eam . tmpo'ing the 'harp diamond s*av s to vtcate a diamond in the right 
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Figure A.7 The Diffuse Diamond configuration. 

refiection of the input beam within the crystal creating a diamond but there is no 
reflection of the secondary beam off of the rear face of the crystal. Absorption of the 
reflected beam by two-beam coupling with the incident beam within the first few 
micrometer* of the crystal edge is probably responsible for the lack of specular 
refiection in this and other similar cases. There is a similar theory which accounts tor 
'he lack of specular refiection of the phase-conjugate be im off of surfaces such as the 
input face of the cry stal. Pepper (Ref. 36| proposes that the lack of a refiection of the 
phase-conjugate beam is due to destructive interference between the incident beam 
reflection, and a beam created by the phase conjugation of the refiection of the original 
phase-conjugate beam. This has been borne out in experiment. However, as Figure 
\ •< .lemcns'ri'es. it certain ancles *he '•'ecular refiection is sufficient!'. m'c:>c • - 

__ . i-.i.T"..!. " us • • ,:'.c..." . 'lie,’. ;' s 

iac\ or a retiection iiave set to be Iuliy determined 

4. Diagonals 

F.gurc A lu illustrates the diagonals The diagonals faded m and cut as Me 
phase-tonpugatc rerlectr.it> oscillated from 12 ' ■ to 21'’« with a period on the order M 











Figure A.3 Example of lack, of reflection of the incident beam 

within the crystal. 


-Is. T: e reflectivity was strongest when the diagonals were brightest When the 
incidence was increased to 06 ° a similar configuration appeared, but with the 
jr.jug.ue reflectivity oscillating from one to 15' •> and specular reflection off of 
s tace lading as the diagonals became brighter. At another angle of incidence 
.enhguraticn ( 62 ° • it was noted that increasing the power of the inpu* beam 
nun :: e a.l oscillations ot the phase-conjugate reflection while elinunat.ni the 
•' fh;s configuration has also been seen without any measurable nha^c- 

Temporal Effects and Frequencs Shifts 
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Figure A.9 Example of multiple reflections within the crystal. 

/. Amplitude Oscillations 

Frequently, amplitude oscillations were observed in the phase-conjugate 
signal of many of the previously discussed examples. These include Figures A.2. A d, 
A.b, A.3. A.9. The oscillations, probably similar to those reported by Valley and 
Dunning in an external resonating cavity jRef ?9j, were most obvious when a mottled 
'"'earn was focused into the sample, but they also appeared in some configurations with 
a.i beam profiles. The observed oscillations often appeared to be chaotic, and much 
work, has been reported recently concerning the observation of chaotic oscillations in 
BaT.O, [Ref. ?d,?5] However, some cf the observed oscillations were sen. near!’, 
periodic ■ periodic w-.thm Id” ■> of some mean period' with periods ranging from r,|| i to 


w:’:. .i ;.used Gaussian beam incident or. the center of the crystal do 1 ' to the c-ax.s 
1 hose . seulations appeared imrueduteiv upon the initiation of phase comugut :. 


I hese were maintained as shown lor a period of a rout one minute and then i u rup - 
shirked to the oscw.ati.ins shown f cure \ 12 Within a lew minutes die ncu 



















Figure A. 10 The Diagonals configuration 

oscillations reappeared Fast oscillations ( —^10 Hz> amplitude modulated b\ slower 
oscillations 1 "*0.5 Hzh and oscillations with faster beats superimposed on them were 
also observed Nowak. has shown that the frequency of these amplitude oscillations :s 
proportional to the mtensite of the incident beam [Ref ?“• ar.d this has been confirmed 
others Ref 

These amplitude osculations of the phase-con;ugate signal are probable a 
result cf frequency shifts of those au\:liarv -"earns that account tor the formation o: tite 
grut.r.g Frequence sh.fts of the aux.iiare beams were proposed b\ I cm berg ar'.d 
3 a otter Ref 4<> to explain observed trequer.ee sweeping in the output :t a ring 
re* •: ,t y : Bain), and in ar. m.iso.ated Ar*’ laser producing phase *o:v..gat.ot. 
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Figure A. 12 Amplitude oscillations which progressed from 
those of Figure A. 11. 

:c"e',3t:on of amplitude and frequency shift of the phase-conjugate beam may have 
-eer. exp!uned by MacDonald and Feinberg when they showed that the four-wave 
mixing signal fails off smoothly as the frequency difference between the reflected and 
incident beam becomes larger than ~ 1 Hz. Although the mechanism for this 
frequency shifting is still under investigation, Lam [Ref. 3] has proposed that the bulk 
pnctovoltaic effect gives rise to an intensity-dependent frequency shift. This implies 
that the frequency is dependent on the amplitude instead of vice versa. Another 
explanation may be that there is consistent movement of the grating when the reflected 
pumpwave forms a loop (see Figure A.2 for example). In this case each circuit around 
the loop -.v-.;: Dcprier-shift the pumpwave and the observed effect may be expected. 

* .ere v>< mm me imer-eu instance vnere me frequency mierenm 
j m ...e - v. o oeair.s a as non-zero ana constant [no photograph available). The 
phasc-comugatc return was frequency-shifted from the input beam, but there was no 
'.leeping That is. there was obvious movement of the interference fringes at a 
mr.star.t ve.ocitv This elTcct resulted from a mottled Gaussian beam incident on the 
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crystal at an angle of 35° to the c-axis. However this :> ar. jr.onta. -• r 

yet to be repeated or explained. 

3. Influence of Returning the transmitted Beam 

This effect is similar to that reported by f einberg Re: 2" i . v 
discovered while trying to stabilize the amplitude osculations of the y 
signal discussed above. The effect was seen using an arrangement ’ - 

nonperiodic amplitude oscillations of the phase-conjugate reflection ranging • ■ / 

to 46%. In this configuration a 100% reflectance mirror was placed behind t c- m- • ■ : 
directing the transmitted beam back through the crystal, but not direcv . 
incident path. All oscillations stopped, leaving the phase-conjugate reflect. or. 

10 %. 

The explanation for this effect is probably found m the conipevi.. • 
between beam patterns. In any instance, the observed beam pattern is the me .sh.^u 
exhibits the maximum gain. When oscillation occurs between two competing w,.:-. 
patterns it is apparent that, for some reason, the gain of one is falling abuse ar.d tr.er. 
below the possible gain of another. The input of a beam of light from any source a /:. 
brightness on the order of the internal beams will cause a change in the gam oi \.e 
system due to some or all of the possible nonlinear effects. In the case of returning the 
transmitted beam, it appears that the gain of one of the competing patterns is 
sufficiently reduced to eliminate it from the possible patterns, 
d. The Free Standing Auxiliary Beam Pattern 

This section presents an occasionally observed result due to a small crystal 
motion. Here the crystal was translated at right angles to the beam path while the 
internal beam pattern was carefully observed for any changes. 

Note the beam configuration seen in Figure A. 13. When the crystal was 
shifted slightly to the left, the auxiliary beam pattern moved with the crystal as shown 
in Figure A. 14 instead of changing to adapt to the new crystal orientation. Note, 
however, that some slight distortion of the pattern did occur. In this case the phase- 
conjugate reflectivity was originally 32% ’with small amplitude oscillations on the order 
of 5%. There was no evidence jf freuuency shifting of the nhase-conmuate \m::: 
Upon translation of the crystal, the phase-conjugate reflectivity diminished to 23% 
with oscillations of the same amplitude as before the translation. 

After translation, any rotational motion of the crystal destroys the auxiliary- 
beam and a new one appears that is similar to the original pattern (i.e. Figure A 13 > 


















Figure A 14 Free standing beams created by transla'ing the 
crystal ''’0.5mm from the configuration of Figure A. 13. 






APPENDIX B 
PSIM LISTING 
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!PSIM_3 

! A SIMULATION OF INTERNAL BEAM RESPONSE IN BARIUM TITANATE 

! IBEAM » SIZE OF INPUT BEAM 

! IBEAMF * ACTUAL NUMBER OF POINTS IN WAVEFRONT AFTER 
! DESIGNATION OF ANGLE PHI 

! ISIZE ■ SIZE OF CRYSTAL 

! ICRYS () » CRYSTAL GRID 

! KEEP() * ARRAY TO KEEP TRACK OF CURRENTLY OCCUPIED SPACES IN THE 
GRID, 2=COL, l=ROW 
! IPHI ■ INPUT ANGLE OF INCIDENCE 
I PHI ■ ANGLE OF INCIDENT BEAM (REAL) 

I PHIC ■ ANGLE OF INCIDENT BEAM AFTER ENTERING CRYSTAL 
! IX a NUMBER OF COLUMN POINTS ACROSS WAVEFRONT 
IY = NUMBER OF ROW POINTS ACROSS WAVEFRONT 
! ADD() * AMOUNT TO ADD TO Y FOR EACH X MOVED 
! IROW = ROW NUMBER 

! ICOL = COLUMN NUMBER 

! ITOP » HIGHEST INDEX OF REFRACTION 

! IBOT = LOWEST INDEX OF REFRACTION 

! GRAD * GRADIENT OF INDEX OF REFRACTION ACROSS THE BEAM 
! CRYS = REAL VALUE OF INDEX OF REFRACTION (BEFORE FIXING) 

! IMDEX = INDEX OF REFRACTION OF LOCATION 

REALX = PLACE FOR KEEPING TRACK OF AMOUNT TO ADD TO COLUMN 
(OR ROW FOR REALY) WHEN MOVING DOWN THE BEAM 
THIS ENSURES THAT ALL FRACTIONS OF A GRID SPACE 
ARE NOT LOST 
REALY = SEE REALX 

REALAD () = SAME AS REALX ABOVE BUT ONE SPACE FOR EACH POINT. 

USED DURING ACTUAL WAVEFRONT MOVEMENT 
IPERMC() = STARTING COLUMN FOR WAVEFRONT POINT 

IMARKB aMARKER TO KNOW WHEN REFLECTION OFF OF BOTTOM HAS OCCURRED 
IMARKS = MARKER TO KNOW WHEN REFLECTION OFF OF SIDE HAS OCCURRED 
ICHECK a INDEX OF TIME PERIODS 
RSTOP a NUMBER OF TIME PERIODS IN SIMULATION 


IDIR » NUMBER OF TIME PERIODS BETWEEN COMPUTATIONS 
OF DIRECTION OF WAVEFRONT 
MDIR = COUNTER FOR CHECKING AGAINST IDIR 
XX = HORIZONTAL DISTANCE BETWEEN END POINTS OF WAVEFRONT 
w = 


: :hnc 


VERTICAL DISTANCE BETWEEN END POINTS OF WAVEFRONT 
= VARIABLE USED TC IHANGE DIRECTION AFTER REFLECTION 
- IRAC; , = NAIF CM -JSED TO MARK THE INCIDENT BEAN PATH 50 THAT 
THE INDICES SET UP BY THE INCIDENT BEAM DOMINATE 
MEAN * MEAN VALUE OF THE INDEX OF REFRACTION 
REF = ACTUAL INDEX OF REFRACTION OF THE PHOTOREFRACTOR BEING 
SIMULATED 

IPLACE = POINT THAT BEAM STARTS IN CRYSTAL 

IGRAPH » MARKER FOR POINTS TO BE PLOTTED (SET FOR 100 POINTS) 
IGRMAX » NUMBER OF ITERATIONS BETWEEN PLOTTING WAVEFRONTS 
CAXIS « FRACTION OF TOTAL INDEX CHANGE DUE TO ANGLE BEAM 
MAKES WITH THEE C-AXIS IN THE CRYSTAL 
CINDEX » REAL VALUE OF INDEX OF REFRACTION USED TO COMPUTE 
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840 
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390 
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910 

920 

930 

940 

950 
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970 
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1100 

1120 
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.150 

. 1"0 

.130 
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! THE CHANGE DUE TO THE ANGLE TO C-AXIS 

! IOUTO ■ MATRIX FOR HOLDING VALUES OF PLOTTED POINTS 

! ICOUT * COUNTER FOR PLACE OF POINT IN OUTPUT FILE 

1 ISKIP = MARKER TO CHECK IF WAVEFRONT IS UNDERGOING REFLECTION 
! REFLCHKB “MARKER TO KNOW WHEN REFLECTION OF WAVEFRONT STARTS/ENDS 
I REFLCHKS * SAME BUT FOR SIDE 

! OP1 * NUMBER OF POINTS IN WAVEFRONT TO LOOK AT FOR RECALCULATION 
1 THIS EQUALS IBEAMF/PIECE 

! ISETMARK = MARKER TO TELL WHEN TO STOP LOOKING FOR REFLECTION 

! PIECE = NUMBER OF PIECES OF WAVEFRONT 

! ADDY() ■ FRACTION OF ONE BOX THAT EACH POINT STARTS OUT AHEAD OF 

! ITS ROW 

1 
I 

i 


INTEGER Ibeam, Ibeamf,Isize,Icrys(500,500) 

INTEGER Keep(2.50).Iphi,lx,Iy,I row,Icol,I top,Ibot,Reflchks(50) 
INTEGER Reflchkb(50), Ikfac,This.That 

INTEGER Imdex,Imaexp(50),Ipermc(50),Imarkb(50),Itime(SO),Iskipb 
INTEGER Imarks(50),idir,Mdir,Ichng,Icout,Iskips,Itemp,Opl,Piece 
INTEGER Itrac(500,500).Mean,Iplace,Igraph.IoulTl,50,101),Isetmark 
REAL Phi,Phic,Add(50),Grad ( Crys,Realx,Realad(50),Rstop,Grmax 
REAL Xx,Yy,Ref,Caxis(50),Cindex,Rcheck,Kfac,Realy,Addy(50) 


Iskipb 

.Piece 


! 


! NOTE THAT IF BEAM SIZE EXCEEDS 50 YOU MUST REDIMENSION 
! 

ON ERROR GOSUB 5430 ! WRITES FILE TO DISK IF ERROR 

PRINTER IS 1 ! SET TERMINAL OUTPUT 

GRAPHICS OFF 

OUTPUT KBD; K; !CLEARS SCREEN 

ON KEY 2 LABEL SAVE &STOP GOSUB 5430 
PRINT TO STOP PROGRAM AND SAVE DATA PRESS f2 


INPUT FACTORS 


INPUT NAME OF DATA FILE TO BE USED ,Name$ 

PRINT DATA FILE: ,NameS 

INPUT SIZE OF CRYSTAL TO BE USED ,Isize 

PRINT CRYSTAL SIZE IS ,Isize, ON A SIDE 

INPUT BEAM SIZE IS ,Ibeam 

PRINT WITH A BEAM OF SIZE , Ibeam 

INPUT PLACE ON CRYSTAL FOR BEAM TO ENTER IS , Iplace 

PRINT ENTERING AT , Iplace 

INPUT HIGH VALUE OF INDEX IS , I top 

PRINT LOW VALUE OF INDEX IS ,Ibot 

PRINT INDICES RANGE FROM ,Itop, TO ,Ibot 

INPUT INITIAL ANGLE IN DEGREES IS (89 DEG MAX) , Iphi 

T1=TIMEDATE 

! 

i 

INITIALIZE "ARIA3LZS 

Acnecx - 3 
Mdir=0 
Ref=2.4 
lcout=0 
lsetmark=0 

Piece=5 ! REM TO CHANGE KK IF YOU CHANGE PIECE 
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1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1350 
1360 
. 37C 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

1960 

1970 

1980 


PRINT NUMBER OF PIECES IN WAVEFRONT , Piece 
1 
! 


» 

i 

! 


PUT FACTORS IN OUTPUT FILE 


Iout(0,0,0)*Isize 
Iout(0,l,0)*Ibeam 
lout(l,0,0)*lplace 
lout (1,1,0)*Itop 
lout(0,0,101)*lDot 
lout(0,1,10l)*Iphi 
lout(1,0,101)“Piece 
lgraph*0 

Rstop«Itop*2 I SET TIME TOTAL PERIODS 

Rstop*Rstop*Isize !USE TWO LINES TO KEEP FROM INTEGER PROBLEMS 
PRINT TOTAL TIME * , Rstop 

Grmax-Rstop/100 ! SET FOR 100 WAVEFRONTS PLOTTED 
GOSUB 5170 !SET UP GRAPHICS DISPLAY 
! 

I TIME BETWEEN NONLINEAR CALCULATIONS 
Idir*5*Itop 
! 

! 

I INITIALIZE ARRAY 
J 

FOR 1*1 TO Ibeam 


>(l,I)-0 
>( 2 , I)*0 


Keep 
Keep( . 

Itime(I)=0 
Reflchks(l)=0 
Reflchkb(l)=0 
NEXT I 
! 

! 

! 

! CONVERT PHI FROM DEGREES TO RADIANS AND FIND ANGLE INSIDE CRYSTAL 
Phi*Iphi 

Phi=Phi*(2*PI/360) 

Phic=(SIN(Phi)/Ref) 

Phic=ASN(Phic) 


CHANGE SIZE OF BEAM DEPENDENT ON THE INPUT ANGLE 


X=Ibeam*COS (Phic} 

Y=Ibeam*SIN(Phic) 

Iy=Y 

Ix=X 

Ibeamf=Ix 

Opl=Ibeamf/Piece IISIZE OF PIECE TO LOOK AT WHEN CALCULATING DIRECTION 

Add(1)=Y/X 

Addl*Add(l) 

IF Iplace>(Isize-Ibeamf) THEN ! MAKE SURE IPLACE IS NOT TOO LARGE 
Iplace =Isize-Ibeamf 

END IF 


! 

! 

! 

Irow=0 

Realy=0 

Icol=Iplace 


FILE WAVEFRONT COORD 
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1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2530 

2590 

2500 

2o23 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 


FOR 1 = 1 TO Ibeamf 
Ketp(2,I)=Icol 
ReaIy=Realy+Add(1) 

Add(l)«Addtl) 

IF Realy>Irow+l THEN Irow*Irow+l 
Addy(I)=Realy-Irow 
Keep(l,I)-Irow 
Icol*Icol+l 
NEXT I 


*****Set UP INDEX OF REFRACTION FOR PATH OF INCIDENT BEAM **** 


SET UP CRYSTAL TO MEAN INDEX 

Mean»ABS((Itop+Ibot)/2) 

FOR 1*1 TO Isize 

FOR J-l TO Isize 

Icrys(I,J)*Mean 
Itrac(I,J)*0 
NEXT J 
NEXT I 
! 

I 

! 

Grad*(Itop-Ibot)/2 !INDEX OF REFRACTION GRADIENT 
FOR 1*1 TO Ibeamf 
Irow=Keep(l,I) 

Icol*Keep(2,I) 

Crys*Grad*SIN(((Ibeamf/2)-(I-l))*(PI/Ibeamf)) 
Crys=Crys+Mean 
Icrys(Irow,Icol)*Crys 
NEXT I 
1 
! 
i 

i SET UP VALUES DOWN THE BEAM PATH 
! 

FOR 1=1 TO Ibeamf 
Icol=Keep(2,I) 

Irov=Keep(l,I) 

Imdex=Ic rys(Irow,Icol) 

Imdexp(I)=Imdex 

Itemp=icol 

Realx=0 

FOR J*Irow TO Isize 

Icrys(J,Itemp)*Imdex 
Itrac(J,Itemp)=l 
Realx=Realx+Aad(I) 

Itemp=Icol-Realx 


uex: 


:jezt 


IF Itemp>Isize THEN Itemp 

I? I , :emD<0 THEN ItemD=0 


=Isize 


SET THE STARTING COLUMNS FOR THE WAVEFRONT AND 
INITIALIZE FLAGS 


FOR 1=1 TO IBEAMF 

Ipe rmc(I)=Keep(2,1) 
Imarkb(I)=0 







2720 

2730 

2740 

2750 

2740 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 

3140 

3150 

3160 

3170 

3180 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3290 

3200 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3390 

3400 

3410 


I»arka(I)«0 
HI XT I 


MOVEMENT *********************************** 

MAKE SURE BEAM IS STILL IN CRYSTAL 

FOR 1*1 TO Ibaamf 

IF Imarkb(I)-l THEN 

IF Kaap(l,I)<2 THEN GOTO 5140 
END IF 

IF Imarks(I)»l THEN 

IF K««p(2,I)>Isize-10 THEN GOTO 5140 
END IF 
NEXT I 
1 

j *******movement ****************************************** 
! 

1 

WHILE Rcheck<Rstop 
1 

FOR I»1 TO Ibeamf 
Icol«Keep(2,I) 

Irow»Keep(l.1) 

Imdex*Icrys(I row,Icol) 

! 

! MARK REFLECTION IF NECESSARY 
1 

IF Icol<=l THEN 

IF Imarks(I)=0 THEN 
Imarks(I)*l 
Realad(I)*l 
END IF 
END IF 

IF Irow =Isize THEN 
Imarkb(I)=l 
END IF 

! 

! 

! MOVE IF IT IS TIME 
! 

IF Itime(I)<Imdex THEN 
Itime(I)=Itime(I)+l 
ELSE 

Itime(I)=0 

Realaa(i)=Realad(I)+Add(I) 

IF Imarks(I)=0 THEN 

Icol=Ipermc(I)-Realad(I) 

IF IcolcO THEN lcol=0 
ELSE 

Icoi=Realad(I) 

END IF 

Keep^.2,10=lcol 
IF Imarkab(I)=0 THEN 
Keep(l,I)=Irow+l 

Keep(l,I)=Irow-l 
END IF 


END IF 
NEXT I 




3420 1 

3430 1 PUT POINTS IN FILE AND DISPLAY GRAPHICS 

3440 l IF IT IS TIME TO DO SO 

3450 1 

3460 I 

3470 Igraph«Igraph+l 
3480 IF Igrapn>Grmax THEN 
3490 Icout»Icout+l 

3500 FOR L»1 TO Ibeamf 

3510 lgraph«0 

3520 lout(0,L,Icout)*Keep(2,L) !X COORD OF POINT 

3530 Iout(l,L,Icout)=Isize-Keep(l,L) !Y COORD OF POINT 

3540 PENUP 

3550 PLOT Iout(0,L,Icout),Iout(l,L,Icout) 

3560 NEXT L 

3570 END IF 
3580 I 
3590 I 
3600 ! 

3610 I 
3620 t 

3630 Rcheck*Rcheck+l 
3640 t 

3650 ! CHECK IF TIME TO RECOMPUTE DIRECTION ****************** 

3660 1 

3670 MDIR=Mdir+l 
3680 IF Mdir>*Idir THEN 
3690 Mdir=0 

3700 PRINT TABXY(1,11), CURRENT TIME PERIOD ■ , Rcheck 

3710 ! 

3720 I 
3730 ! 

3740 ! 

3750 1 

3760 I 
3770 I 

3780 FOR K-0 TO Piece-1 
3790 ! 

3800 ! **************** SET FOR 5 PIECES ***************** 

3810 ! 

3820 IF K*0 THEN Kk=0 

3830 IF K*1 THEN Kk=4 

3840 IF K-2 THEN Kk=l 

3350 IF K»3 THEN Kk=3 

3860 IF K»4 THEN Kk=>2 

3870 !IF K=5 THEN Kk=7 

3880 I IF K=6 THEN Kk*3 

3890 !IF K=7 THEN Kk=6 

3900 I IF K=8 THEN Kk»4 

3910 !IF K=9 THEN Kk*5 

3920 This=Kk*Opl+l 

3930 That=(Kk+l)*Opl 

3940 IF Kk*Piece-l THEN That*Ibeamf 

3950 ! 

3960 ! 

3970 ! FIRST CHECK TO MAKE SURE THE FRONT IS NOT PARTIALLY REFLECTED 

3930 ! 

3990 : 

4000 -SXipb=0 

4010 I jet.Tiark=0 

4020 IF Reflchkb(K)=0 THEN 

4030 FOR M=This TO That 

4040 IF Imarkb(M)*l THEN 

4050 Isetmark=l 

4060 FOR N=This TO That 

4070 IF Imarkb(N)*0 THEN Iskipb*l 

4080 NEXT N 

4090 END IF 

4100 NEXT M 

4110 END IF 











4120 
4130 
4140 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 
4310 
4320 
4330 
4340 
4350 
4360 
4370 
4380 
4390 
4400 
4410 
4420 
4430 
4440 
4450 
4460 
4470 
4480 
4490 
4500 
4510 
4520 
4530 
4540 
4550 
4560 
4570 
4530 
4590 
4600 
4610 
4620 
4630 
4640 
4650 
4660 
4670 
4630 
-o9 : 

T J * 

"t > 1 <J 

4720 

4730 

4740 

4750 

4760 

4770 

4780 

4790 

4800 

4810 


IF Isetmark*l AND Iskipb-0 THEN Reflchkb(K)*l 

lskips*0 

lsetmark*0 

IF Reflchks(K)*0 THEN 
FOR M*This TO That 
IF Imarks7n)«l THEN 
FOR N*This TO That 
Isetmark*l 

IF Imarks(N)*0 THEN I skips*1 
NEXT N 
END IF 
NEXT M 


IF Isetmark*! AND lskips*0 THEN Reflchks(K)*!!SKIP THIS IF NOT NECESSARY 


END IF 


IF lskips*0 AND lskipb*0 THEN I IF WAVEFRONT IS NOT PARTIALLY REFLECTS! 


SLOPE OF PIECE OF WAVEFRONT IS YY/XX 


y*(Keep(l,This)+Addy(This))-(Keep(l,That)+Addy(That)) 
1=1 time(This)/Icrys(Keep(l.This),Keep(2,This); 


Y2=I time (That J/lcrys (Keep^ 1That))Keep(2 ^That)) 
Y^*Yy+Yl-Y2 


I 


Imarks(This)*0 THEN 

Xx*(Ipermc(That)-Realad(That))-(Ipermc(This)-Realad(This)) 

ELSE 

Xx*Realad(That)-Realad(This) 

END IF 

IF Xx>l THEN 
FOR G-This TO That 

Add(G)*ABS(Yy/Xx) !FIGURE NEW DIRECTION 
Caxis(G)*Xx/Opl 

IF Kk-Piece-1 THEN Caxis(G)*Xx/(Ibeamf-(K*Opl)) 

NEXT G 

PRINT Add(G-l) 

END IF 


I 


! 


FOR I*This TO That 
Icol*Keep(2,I) 
Irov«Keep(l,I) 
Imdex*Imdexp(I) 


; teirp* .. 
.Aea.k=0 


Oi. 


I CHECK FOR REFLECTION 


! 


IF Imarkb(I)*0 THEN 

FOR J*Irow TO Isiae 

Icrys(J,Itemp)*Imdex 
Realx*Realx+Add(I) 

I temp*Icol-Realx 
IF lcemp<0 THEN ltemp*0 




a 


a 




IF Itemp>Isize THEN Itemp*Isize 
Itrac(J,Itemp)*l 
NEXT J 


!INCIDENT BEAM DOMINATES 
FOR J»0 TO I row 

Ichng*Irow-J 

IF Itrac(Ichang,ltemp)*0 THEN 
Icrys(Ichng,Itemp)=Imdex 
ELSE 

Icrys(Ichng,Itemp)*Icrys(Ichng,Itump) 
END IF 
J 

Realx=Realx-Add(I) 

Itemp*Icol-Realx 
IF ItempcO THEN ltemp«0 
IF Itemp>Isize THEN Itemp»Isize 
NEXT J 


END IF 
NEXT I 
END IF 
NEXT K 
END IF 


I CHECK FOR END OF RUN ******************** 

END WHILE 

T2*TIMEDATE 

PRINT TOTAL RUN TIME * ,T2-T1 
GOSUB 5360 !WRITE FILE TO DISK 

STOP 

! SUBROUTINE TO INITIALIZE GRAPHICS 

GRAPHICS ON 

GCLEAR 

SHOW 0, Isize,0,Isize 
PLOT 0.0,-2 
FOR Q-6 TO Isize 
PLOT 0,0,-I 
NEXT Q 

FOR Z»0 TO Isize 
PLOT Isize,Z,-I 
NEXT Z 

FOR Q»Isize TO 0 STEP -I 
PLOT Q, Isize,-1 
NEXT Q 

FOR Z*Isize TO 0 STEP -1 

PLOT O.Z.-I 

NEXT Z 

PLOT 0,0,-2 

RETURN 

! SUBROUTINE TO OUTPUT FILE TO DISK 
MASS STORAGE IS 700,0 
CREATE BDAT Name? 10404,2 
ASSIGN Par.'. I TO 'Jarre5 


! SUBROUTINE TO OUTPUT FILE IF ERROR OCCURS 

MASS STORAGE IS s.700.0 

CREATE BDAT NameS,10404,2 

ASSIGN Path 1 TO NameS 

OUTPUT Path l;Iout(*) 

ASSIGN Path 1 TO * 

PRINT ERROR! PROGRAM TERMINATED 
T2-TIMEDATE 

PRINT TOTAL RUN TIME ■ .T2-TI 
STOP 


1 
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